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1.0   INTRODUCTION 

1. 1  Problem  Formulation  and  Engineering  Significance 

Mankind  has  been  using  organic  fuels  for  industrial  and  domestic 
purposes  for  many  years.   The  fluctuating  price  of  petroleum  and  the 
numerous  hazards  involved  in  the  inefficient  combustion  of  fossil  fuels 
have  motivated  the  diversification  of  fuel  sources  as  well  as  the  U.S. 
Government  and  private  companies  interest  for  new  research  involving  the 
causes  of  inefficient  combustion  in  some  of  the  most  widely  used 
hydrocarbons.   Extensive  and  recent  research  on  fuels  such  as  benzene, 
acetylene*  cyclohexane  etc.  has  demonstrated  that  soot  formation  is  a 
problem  in  the  combustion  of  hydrocarbons.   This  problem  results  in: 
1)  environmental  hazards  from  emission  of  particulates  and  associated 
compounds,  2)  engineering  difficulties  in  the  materials  due  to  the 
thermal  radiation  from  soot  particles,  and  3)  mechanical  or  heat 
transfer  hindrance  caused  by  layers  of  soot  deposited  on  critical 
surfaces  within  a  mechanical  device. 

The  environmental  hazards  that  soot  represent  are  very  well 
known.   Early  papers  in  medicine  such  as  the  one  published  in  England  by 
Sir  Percival  Pott  in  1775  discussed  the  cancer  of  the  scrotum  as  an 
occupational  disease  of  chimney  sweeps.     These  chimneys  contained 
deposits  of  soot  formed  from  the  burning  of  coal.  More  recent  studies 
have  confirmed  that  solid  particles  from  coal  and  oil-fired  power  plants 
contain  a  large  number  of  carcinogenic  compounds.   The  Identification  of 
carcinogenic  compounds  in  soot  has  yielded  a  large  list  of  benzene 
extractable  polycyclic  aromatic  hydrocarbons  (PAH) . 


A  second  problem  that  the  emission  of  soot  represents  is  the 

engineering  difficulties  caused  by  the  thermal  radiation  of  soot 

(2) 
particles.   An  example  of  this  is  discussed  by  Vaughn    who  reported 

that  soot  particles  accumulate  in  the  blade  section  of  an  advanced 

design  high  temperature  gas  turbine  and  cause  detrimental  effects  by  the 

increase  in  temperature  due  to  the  thermal  radiation  from  soot 

particles. 

Finally,  another  problem  where  soot  formation  can  be  held 
responsible  is  the  mechanical  or  heat  transfer  inefficiency  of 
mechanical  devices  by  soot  proliferation  on  critical  surfaces. 

All  of  these  undesirable  consequences  of  soot  formation  have 
produced  a  special  interest  in  the  pyrolysis  of  aromatic  hydrocarbons 
for  the  present  research  work,  where  toluene  was  chosen  to  increase  our 
understanding  of  the  initial  steps  involved  in  the  production  of  soot. 

1.2  Thermal  Decomposition  of  Toluene 

The  present  section  gives  a  brief  description  of  the  experimental 
work  performed  on  toluene  pyrolysis.   The  discussion  format  starts  with 
a  background  information  about  toluene;  then,  special  attention  is  given 
to  the  different  experimental  studies  on  toluene  pyrolysis  reported  in 
the  literature. 

The  discovery  of  toluene  goes  back  to  the  mid  1830's  when  Pelletier 

(3) 
and  Walter    were  able  to  obtain  this  product  by  heating  rosin  oil,  and 

(4) 
later  on,  Deville    reported  his  results  on  the  distillation  of  the 

Balsam  of  Tolu,  a  natural  resin  shipped  to  Europe  from  the  seaport  of 

Tolu  in  Colombia,  South  America. 
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Toluene  was  officially  recognized  as  a  new  hydrocarbon  later  in 
1847,  when  Charles  Mansfield  separated  this  product  and  some  others  from 
coal  tar.     Despite  the  fact  that  toluene  was  quickly  recognized  as  a 
good  solvent,  significant  commercial  use  did  not  occur  until  it  was 
realized  that  the  nitration  of  toluene  would  produce  T.N.T. 
(Tri-nitro-toluene) ,  a  very  powerful  explosive.   It  is  also  a  solvent 
for  paint  and  ink  and  used  in  printing  newspapers  and  magazines.   One  of 
its  most  important  applications  came  with  the  invention  of  the  internal 
combustion  engine  where  it  has  been  used  with  other  aromatic 
hydrocarbons  in  gasoline  to  suppress  engine  knock,  thereby  increasing 
the  power  and  providing  a  smoother  running  engine. 

Toluene  pyrolysis  has  been  a  subject  of  interest  in  combustion 
studies.   These  studies  can  be  divided  into  three  different  areas  due  to 
the  treatment  that  experimentalists  have  given  to  its  pyrolysis:  1)  soot 
formation  studies,  2)  kinetic  studies,  and  3)  product  yield  studies. 

1.2.1  Soot  Formation  Studies 

Studies  on  soot  formation  from  toluene  pyrolysis  are  relatively  a 
new  subject  that  experimentalists  have  been  concentrating  on  due  to  the 
similarity  in  the  soot  yields  with  other  aromatic  hydrocarbons  (such  as 
benzene).   Three  theories  on  soot  formation  from  pyrolyzed  hydrocarbons 


have  been  proposed  from  studies  as  early  as  1953.   '    These  theories 

(2) 
are  also  discussed  by  Vaughn    and  they  are:  1)  The  acetylene  theory, 

which  proposes  that  a  decomposed  hydrocarbon  forms  unsaturated  fragments 

(such  as  acetylene)  that  polymerize  and  produce  larger  species  that 

eventually  form  soot.   2)  The  C_  polymerization  theory,  which  proposes 

that  the  pyrolyzed  hydrocarbon  forms  C  intermediate  fragments  that 

represent  nucleation  sites  for  polymerizations  with  other  species  or  the 

C.  themselves  polymerize  to  form  soot.   And  finally,  3)  the 

polymerization  theory  that  proposes  a  series  of  polymerization  reactions 

to  form  a  high  molecular  weight  polymer  that  dehydrogenates  to  form  soot 

particles. 

A  general  mechanism  that  most  authors  have  adopted  has  been  the  one 

(8) 
proposed  by  Graham,  et  al.    working  with  a  shock  tube  in  the  high 

temperature  region  1600-2300  K  (reaction  time  2.5  msec).   Experiments 

with  benzene,  toluene,  ethylbenzene  and  indene,  highly  diluted  with 

argon  led  to  the  conclusion  that  soot  is  formed  from  aromatic 

hydrocarbons  in  two  different  pathways:  1)  a  direct  fast  route,  with 

condensation  reactions  and  almost  instantaneous  soot  formation  at 

temperatures  below  1750  K,  and  2)  an  indirect  slow  route  from  the 

hydrocarbon  fragments  produced  at  high  temperature.   This  mechanism  is 

shown  below  and  it  condenses  the  three  proposed  theories  for  soot 

formation  outlined  earlier. 


CH.,     Condensation  ^\  >^\.       fast 


Soot 


reactions    N^^-v^-^  (direct  route) 


CH, 


Fragmentation       CH 

*    C2K 

reactions       C,H 
3  x 


slow 


■*   Soot 

(indirect  route) 


Although  Graham,  et  al.  considered  that  temperature  was  the 
determinant  parameter  for  the  soot  formation,  other  studies  have  proved 

that  reaction  time  is  also  a  parameter  of  important  consideration  in  the 

(9) 
proposal  of  soot  formation  routes.   For  instance,  Fields  and  Meyer son 

pyrolyzed  toluene  and  toluene-ct,d,  at  low  temperature  (873  K)  and  large 

contact  time  (6  sec) .   They  proposed  a  bimolecular  complex  via  such 

reactions  as 


2<V%CHS— 


where  toluene  can  form  bibenzyl  via  dlhydrogenation  or  can  collapse  to  a 
benzyl-or  tolyl-methyl  cyclohexadiene  which  loses  hydrogen  to  yield 
methyldiphenylmethane  or  dimethylbiphenyl.   Similarly,  Kershaw     in 
1978  working  in  a  rotating  autoclave  at  low  temperatures  (623  K  -  723  K) 
and  long  reaction  times  (1-3  hours)  found  a  high  concentration  of 
bibenzyl,  diphenylmethane  and  methylbiphenyl  together  with  three  isomers 
of  bibenzyl,  which  are  also  soot  precursors.   In  the  present  work,  it 


will  be  seen  later  that  these  high  molecular  products  were  observed  in 
lower  concentration  due  to  the  low  reaction  time. 

In  1981,  Evans  and  Williams     working  with  shock  tube  and 
mixtures  of  toluene  and  argon  (C:Ar  ratio  0.01:1)  reported  soot  yield 
curves  for  the  temperature  range  715-2159  K  where  the  soot  yield  was 
zero  at  around  1000  K,  increased  to  60%  at  1200  K  and  then  decreased  to 
30%  at  around  2200  K  (see  Fig.  1).   The  reaction  time  in  this  study  was 
not  reported,  but  from  results  in  the  present  study  it  is  presumed  to  be 
much  greater  than  2  msec.   Evans  and  Williams  agreed  with  the  mechanism 
for  soot  formation  proposed  by  Homann  and  Wagner     that  consisted  of 
the  step-wise  addition  of  C.   radicals  of  acetylene  to  form 
polyacetylenes.   The  polyacetylene  radicals  combine  with  other 

polyacetylenes  to  produce  large  branched  radicals.   The  acetylene 

(12) 
radicals  are  formed  from  the  pyrolysis  of  acetylene  and  ethylene. 

Studies  on  soot  formation  from  toluene  by  Wang,  Ma tula  and 

(13) 
Farmer     using  a  shock  tube  (1981)  produced  some  new  and  interesting 

results.   They  concluded  that  condensation,  polymerization  and 

dehydrogenation  are  the  most  important  soot  formation  steps  not  only  for 

toluene  but  also  for  other  fuels,  and  fragmentation  seemed  to  have  a 

less  important  role.   From  their  experimental  results  they  also 

concluded  that  at  temperature  lower  than  1800  K  an  increase  in  toluene 

concentration  increased  the  soot  yields  due  to  an  acceleration  of  the 

condensation  process.   Hydrogen  and  oxygen  were  introduced  to  the 

system,  and  it  was  observed  that  an  increase  in  oxygen  or  hydrogen 

concentration  in  the  toluene/argon  mixture  would  decrease  the  soot 

yields.   At  low  temperature,  the  suppression  effect  of  oxygen  (or 


Fig.  1.   Variation  of  relative  Soot  Concentration  with 
Temperature  during  the  Pyrolysis  of  Toluene. 
(Evans  and  Williams).    ' 
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hydrogen)  in  soot  formation  was  smaller  than  at  higher  temperature,  this 
was  explained  by  Wang,  et  al.  arguing  that  the  oxygen  (or  hydrogen) 
molecule  attacked  the  methyl  group  at  low  temperature  allowing  for 
condensation  of  phenyl  groups  to  form  soot.   However  at  high 
temperature,  oxygen  would  attack  the  phenyl  radical  suppressing  the  soot 
formation. 

The  latest  studies  of  soot  formation  from  toluene  were  performed  by 
Frenklach  et  al.   4*15'  in  1983  and  1984.   In  these  studies,  the 
combustion  of  toluene  and  oxygen  mixtures  in  argon  was  observed  behind 
reflected  shock  waves  in  the  temperature  range  1496-2391  K  and  pressure 

range  1.85-3.04  atm.   The  mol  fraction  of  toluene  in  argon  was  changed 

-3  -2 

from  3.11  x  10   to  1.75  x  10   and  the  toluene/oxygen  ratio  was 

maintained  at  one.   The  results  ratified  earlier  studies  by  Wang,  et  al. 

where  they  found  that  an  increase  In  toluene  concentration  at  low 

temperature  increased  soot  yields.   However,  Frenklach  et  al.  explains 

this  behavior  by  stating  that  "addition  of  oxygen  causes  reactions  to 

occur  that  compete  with  pressure-dependent  fragmentation  of  the  aromatic 

ring"  which  is  similar  to  the  physical  explanation  given  by  Wang  et  al. 

Also,  Frenklach  et  al.  presented  an  empirical  modeling  of  soot 

formation     for  given  reaction  times  at  various  temperatures, 

pressures  and  initial  reactant  concentrations  which  could  be  used  for 

different  fuels.   The  model  included  fragmentation  and  polymerization 

processes,  and  soot  yields  from  toluene/argon  mixtures  given  the 

following  conditions : 


Temperature  range:  1500-2400  K 

Pressure  (P.):  2.05  to  2.98  atm 

Time:  O.5.-  2.0  msec. 

Toluene  mol  fraction  in  Argon:  10   to  10 


The  results  for  three  different  cases:  pressure  dependence,  variation  of 
reaction  time  and  concentration  dependence,  were  compared  to 
experimental  results  and  the  predictive  power  of  the  model  was  in  good 
agreement  with  the  experimental  data.   Figure  2  shows  the  experimental 
data  and  predicted  results  from  the  model. 

1.2.2  Kinetic  Studies 

The  earliest  kinetic  studies  in  toluene  pyrolysis  performed  by 

(39) 
Tilicheev  (1939)  and  M.  Swarc  (1947)     reported  values  for  the  bond 

energy  of  the  C-H  bond  of  293.1  kJ  and  324.5  kj  respectively.   Also, 

Swarc  found  that  the  reaction  is  of  first  order  by  changing  the  contact 

time  and  the  toluene  pressure  in  the  system.   The  kinetic  rate  constant 

that  he  reported  for  the  initial  step  reaction  in  toluene  decomposition 

C7Hg  ►  C?H7  +  H-  (1.1) 

was  ^  =  2.0  x  1013  exp(-324.5  (kJ)/RT)  sec-1. 

A  higher  value  for  the  C-H  bond  energy  was  reported  by  H.R. 
Anderson,  et  al.     as  374.7  kJ  at  room  temperature  working  on  the 
bromination  of  toluene.   However,  the  discrepancy  with  early  results 
from  Swarc  was  not  discussed  and  was  only  left  to  future  studies  on  this 
matter. 
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Fig.  2. 


Experimental  Data  and  Predicted  Values  from  Empirical  £Iodel 
for  Soot  Formation  from  Toluene  by  Frenklach,  et  al. 

(a)  Constant  reaction  time,  varying  concentration  and 
Pressure. 

(b)  Constant  concentration,  varying  reaction  time. 

(c)  Constant  reaction  time,  varying  concentration. 
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Kinetic  studies  by  Blades,  Blades  and  Steacie     were  not 
consistent  with  previous  observations  by  Swarc.   This  newer  study  found 
that  the  first  order  rate  varied  with  pressure,  contact  time,  and  the 
experimental  surface  reactor  vessel.   An  activation  energy  of  376.8  kj 
and  a  frequency  factor  of  10   "  were  reported,  so  postulation  of  the 
equation 


c7h8 k:6h5  ♦  CH3-  (1.2) 

with  reaction  (1.1)  was  assumed  to  be  possible.   This  postulation  and 
the  deviation  from  first  order  behavior  made  it  "impossible  to  assign 

the  observed  activation  energy  to  a  specific  step  in  the  process." 

(18) 
Experimental  results  from  M.  Takahasi     found  the  log  K.  vs  1/T 

plots  to  be  concave  (!L  is  the  first  order  reaction  rate  constant) , 

therefore,  the  temperature  range  (1011-1206  K)  was  divided  into  two 

regions.   It  was  observed  that  the  reaction  rate  constants  for  the  two 

temperature  regions  varied  as  contact  time,  volume  and  surface  area  of 

the  reaction  vessel  were  changed.   This  behavior  was  presumed  to  be  due 

to  the  suppression  of  any  unfavorable  secondary  reaction;  however,  there 

was  no  experimental  evidence  that  supported  this  theory.   Takahasi  then 

discussed  three  possible  causes  of  the  concave  behavior  in  the  log  K. 

vs.  1/T  plots:   1)  experimental  errors  of  the  small  amount  of  products 

present  at  lower  reaction  temperature,  2)  correction  used  for  reaction 

rate  constant  using  a  more  elaborate  Arrhenius  form 

K  =  A  Tn  exp(-E/RT)  (1.3) 

where  A  is  preexponential  term,  T  is  temperature,  n  is  a  temperature 
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correction  factor  assumed  to  be  0.5,  E  is  energy  of  activation  and  R  is 
the  universal  gas  constant,  and  3)  complications  arising  from  competing 

reactions  with  different  activation  energies. 

(19) 
In  1961,  Price     found  that  the  first-order  rate  constants 

decreased  approximately  9%  when  the  contact  time  was  reduced  from  1.0  to 

0.41  seconds.   The  reaction  rate  constant  for  reaction  (1.1)  was  found 

to  be  K  =  1014-8  exp(-347.5  (kj) /RT)  sec"1. 

In  1971,  Peacock,  Brooks  and  Cummins     working  with  a  static 

system  reported  a  value  of  K  =  10      exp(-372.2  (kJ)/RT)  sec"1  for 

reaction  (1.1)  and  an  estimate  for  the  reaction 


C?H7  +  H2  ►  C?H8  +  H-  (1.4) 

9  45  -1 

of  K,  =■  10  '   exp(-60.7  (kJ) /RT)  sec   was  established.   Energies  of 

activation  for  the  formation  of  hydrogen  and  methane  were  also  reported 

as  309.8  kj/mol  and  374.3  kj/mol  respectively.   The  addition  of  hydrogen 

to  the  system  accelerated  the  reaction  process. 

More  recent  kinetic  studies  in  toluene  pyrolysis  performed  in  the 

early  to  mid  1980 's  include  experimental  work  performed  by  Troe, 

(21)  (22)  C23) 

et.  al.,v    Banerjee,  et  al.,     and  Skinner  et.al.    '      Troe  working 

with  a  shock  tube  found  a  limiting  high  pressure  rate  constant 

K1  =  8  x  10   exp(-380.0  (kJ)/RT)  sec"  for  the  temperature  range 

1100-1800  K.   Banerjee  et  al.  reported  a  value  of  K.  =  108"567 

exp(-186.5  (kJ)/RT)  for  the  range  873-1000  K  which  he  recognized  to  have 

a  lower  activation  energy  due  to  "the  different  experimental  conditions 

and  calculation  methods  used".   Finally,  Skinner  et  al.  in  1984  using  a 

shock  tube  calculated  a  limiting  high  pressure  value  of  K.  ■  2.7  x  10 
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exp  (-347.5  (kJ) /RT)  sec   for  reaction  (1.1)  in  the  temperature  range 
1000-1800  K.   This  last  paper  includes  reported  rate  constant  values  for 
other  reactions  occurring  in  toluene  and  deuterated  toluene  pyrolysis. 

1.2.3  Product  Yield  Studies 


This  last  section  in  reported  toluene  pyrolysis  studies  gives  a 

brief  discussion  on  the  product  yields  observed  in  previous  works. 

(24) 
In  1867,  Berthelot     found  hydrogen,  methane,  acetylene,  benzene, 

and  bibenzyl  and  other  aromatic  hydrocarbons  with  condensed  rings. 

Later  in  the  late  1930's  and  early  1940' s,  Tilicheev  and  Swarc  observed 

bibenzyl  as  the  main  product  yield,  and  immediately  reaction  mechanisms 

were  postulated  (see  Results  and  Discussion).   Steacie  et  al.     found 

additional  compounds  such  as  styrene,  dimethyl  biphenyls ,  and  anthracene 

which  were  also  reported  by  Takahasi.      In  1961,  Price     observed 

small  amounts  of  ethylene  and  ethane.   More  recent  studies  performed  by 

R.  SmitlT  '2  '  and  J.R.  Kershaw^10'  report  the  yield  of  high  and  low 

molecular  weight  compounds.   These  compounds  are  shown  in  Tables  1,  2 

and  8.   From  the  studies  performed  by  Smith,  it  is  concluded  for  the 

temperature  range  (1173  -  2073  K)  that  the  highest  yield  product  is 

acetylene  whose  yield  increases  with  temperature,  indicating  the 

disappearance  of  all  compounds  ultimately  to  form  acetylene.   Other  high 

yield  compounds  were  methane  (whose  yield  reaches  a  maximum  of  24%  at 

1573  K),  benzene  (11%  at  1673  K) ,  diacetylene  (6.3%  at  1673  K) ,  and 

ethylene  (5.2%  at  1673  K) . 

The  present  study  on  toluene  decomposition  was  based  on  results 

from  all  these  previous  studies.   The  purpose  of  this  research  was  to 
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investigate  the  thermal  decomposition  of  toluene  in  vapor  phase  in  the 
temperature  range  800  K  to  1700  K  using  a  single  pulse  shock  tube.   The 
yields  of  vapor  phase  and  solid  products  were  measured  and  a  mechanism 
for  the  production  of  the  major  products  was  proposed  using  available 
kinetic  data  from  the  literature. 
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Table  1.   Heavy  Species  Found  . 
by  R.  Smith  (1979) (      ' 


Product 

C8H8 

Styrene 

C8H10 

Ethylbenzene 
p-Xylene 

C10H8 

Naphthalene 

C12H10 

Biphenyl 

C13H10 

Fluorene 

C13H12 

Methyl  biphenyl 

CUH10 

Phenanthrene 

C14H10 

Anthracene 

C14H12 

Stilbene 

C14H14 

Bibenzyl 
Benzyl  toluene 
Dimethyl  diphenyl 

C16H10 

Fluoranthrene 
Pyrene 

C17H12 

Benzof luorene 

C18H12 

Chrysene 

1 , 2-Benz anthracene 

C19H14 

Methylchrysene 

C20H12 

Perylene 

Benzof luoranthene 

Benzopyrene 
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Table  2.   Heavy  Species  Found  . 
By  Kershaw  (1978) X      ' 


Product 

2-Methylbiphenyl 

Diphenylme thane 

3-Methylbiphenyl 

4-Methylbiphenyl 
Bibenzyl 

4-Me thy ldiphenylme thane 
3,3' -Dlmethylbiphenyl 

4,4' -Dlmethylbiphenyl 
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2.0   EXPERIMENTAL  METHODS  AND  TECHNIQUES 

2. 1   Introduction 

The  pyrolysis  of  toluene  has  been  studied  throughout  the  years 

using  different  techniques.   The  first  reported  technique  goes  back  to 

(24) 
1867  when  Berthelot     passed  toluene  through  a  glowing  red  porcelain 

tube  in  which  he  found  similar  products  to  those  found  in  the  present 

work.   Later  on,  Tilicheev  in  the  Soviet  Union  (1939)  and  Swarc  in 

England  (1947)  continued  experimental  work  on  toluene  pyrolysis  using  a 

flow  reactor  which  was  a  very  popular  technique  used  in  the  early 

studies  of  the  pyrolysis  of  toluene  and  other 

u   j     u    (22,19,18,27,28) 

hydrocarbons.  Other  techniques  used  in  toluene 

(20)  (29) 

pyrolysis  studies  were:  static  systems,     laser  heating,     and  shock 

tube_(23,30,31,32,33) 

Following  the  work  performed  by  Tilicheev  and  Swarc  using  flow 
reactors  (and  later  Takahasi,     Blades,  Blades,  Steacie     and 
Price    ),  Brooks,  Cummins  and  Peacock     in  1971  decided  to  use  a 
static  system  because  they  argued  that  In  flow  systems  the  mixture 
composition  could  not  be  varied  easily  and  the  reaction  time  was 
difficult  to  measure  accurately.   Flow  reactors  were  widely  used  in 
early  investigations  due  to  the  fact  that  the  reactant  and  produced 
gases  were  easy  to  handle  and  relatively  large  amounts  of  gas  products 
could  be  collected  for  analysis.   However,  the  availability  of  modern 
gas-chromatographic  techniques  makes  these  advantages  obsolete. 
Flow  reactor  techniques  in  the  form  of  Knudsen  cells  and  rotating 
autoclaves  were  later  used  by  R.  Smith (   *   '  and  J.  R.  Kershaw^10'  but 
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Che  need  for  a  system  that  could  handle  high  temperature  pyrolysis  (flow 

(34) 
experiments  cannot  exceed  the  softening  point  of  quart    )  and  small 

adjustable  reaction  times  in  the  absence  of  oxygen  made  possible  the  use 

of  the  shock  tube. 

Shock  tube  experiments  have  been  conducted  successfully  in  the 
recent  years  in  toluene  pyrolysis.   These  studies  have  been  basically 
intended  to  observe  the  behavior  of  products  and  radical  species  in 
short  reaction  times  with  fast  cooling  periods,  and  also  for  the 
postulation  of  reaction  mechanisms  and  kinetic  parameters  when  soot 
particulates  are  formed  or  are  in  the  process  of  formation. 

In  the  present  study  a  Single  Pulse  Shock  Tube  (SPST)  was  used  in 
which  the  concentration,  reaction  time,  and  reaction  temperature  were 
prescribed  parameters  and  final  species  concentrations,  overall  reaction 
rate,  and  kinetic  mechanism  for  vapor  phase  toluene  pyrolysis  were  to  be 
determined. 

2.2  Materials 


All  the  experiments  in  the  present  study  were  performed  using 
Toluene  (Grade  99.9)  as  a  reactant,  provided  by  Aldrich  Chemical 
Company.   The  purity  of  toluene  was  checked  by  gas  chromatographic 
analysis  and  there  were  some  unidentified  C,  and  C   impurities  that 
accounted  for  0.5%.   This  toluene  was  used  without  further  purification. 
Carbon  Disulfide  (CS.-Purif ied  Grade,  Fisher  Scientific)  was  used  for 
the  extraction  of  solid  products  from  a  liner  that  fit  into  the  test 
section. 

Accessory  equipment  used  for  the  shock  tube  included  two  mechanical 

pumps  (Duo  Seal,  Model  1402)  and  a  water  cooled  diffusion  pump  (Veeco, 
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Model  EP-2W)  that  evacuated  the  shock  tube  to  1.32  x  10   atm.   Also  a 
digital  oscilloscope  (Tektronix,  model  76234)  connected  to  a  timer 
(Fluke,  model  720012)  with  two  pressure  transducers  in  the  test  section 
to  record  the  pressure  history  and  the  time  for  the  shock  to  pass 
through  the  test  section.   The  signal  from  the  transducers  was  amplified 
by  two  linear  amplifiers  (Ortec,  model  410)  and  an  Ortec  (model  456)  HV 
power  supply  was  also  attached  to  the  system.   Pictures  taken  on  the 
oscilloscope  were  recorded  by  a  Tektronix  (model  C5B)  camera. 

2.3  Analytical  Instruments 

The  analysis  of  all  gaseous  and  solid  (or  liquid)  samples  was 
performed  by  three  gas-chromatographs (GC) :  Carle  (model  311),  Perkin 
Elmer  (model  Sigma  3B) ,  Tracor  (model  560) .   These  analyzers  were 
equipped  each  one  with  flame  ionization  detectors  and  the  Carle  GC  had 
the  option  of  a  thermal  conductivity  detector.   All  three  analyzer 
conditions  are  shown  in  Table  3.   These  analyzers  were  all  connected  to 
a  computer  (Perkin-Elmer  Sigma  15)  that  recorded  the  data. 
Identification  of  the  species  present  in  the  compounds  was  made  by 
injecting  known  samples  and  comparing  retention  times.   Quantitative 
results  to  yield  the  concentrations  of  species  in  the  sample  was  made  by 
comparing  the  sample  areas  from  the  analyzers  with  the  areas  obtained 
from  calibration  curves  of  known  species.   This  is  explained  in  the 
following  equation: 


C  MW 

Cs=ATxAsXMir>  «.l) 

C  S 
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where    C  =  concentration  of  sample  species, 

C  =  known  concentration  of  calibrating  species, 

A  =  area  from  GC  in  sample, 

A  ■  area  of  calibrating  species, 

MW  =  molecular  weight  of  calibrating  species, 

MW  ■  molecular  weight  of  sample  species. 

2.4  KSU  Single  Pulse  Shock  Tube  (KSU-SPST) 
2.4.1.   Introduction 

The  shock  tube  is  a  device  in  which  a  plane  shock  wave  is  produced 
by  the  sudden  bursting  of  a  diaphragm  that  separates  two  pressure 
regions:  region  1  with  a  gas  at  low  pressure  (from  0.25  to  0.66  atm)  and 
region  4  with  a  gas  at  high  pressure  (around  19  atm) .   After  bursting 
the  diaphragm  a  compression  wave  is  formed  in  the  low  pressure  gas. 
This  compression  wave  is  accompanied  by  a  sudden  temperature  and 
pressure  pulse  (regions  2,  3)  which  heats  up  the  gaseous  sample  at  the 
end  of  the  low  pressure  region  (test  section).  Simultaneously  and  in 
the  opposite  direction,  an  expansion  or  rarefaction  wave  moves  into  the 
high  pressure  gas  at  the  speed  of  sound  of  this  gas.   The  reflection  of 
the  shock  at  the  end  wall  produces  a  further  rise  in  temperature  and 
pressure  in  the  low  pressure  gas.   However,  the  arrival  of  the  reflected 
expansion  wave  into  this  region  which  will  be  denoted  with  the  subscript 
5,  produces  a  fast  and  uniform  decrease  in  temperature  and  pressure 
which  will  be  explained  later.   See  Fig.  3  for  better  understanding  of 
this  process.   This  whole  process  occurs  in  a  few  milliseconds  and  after 
each  shock  wave  experiment  is  performed,  all  samples  are  collected  to  be 
analyzed  by  gas  and  liquid  chromatography. 
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2.4.2  KSU-SPST  Design  and  Operation 


The  Kansas  State  University  shock  tube  was  first  built  and  tested 

5!. 
(2) 


(35) 
by  Seeker  in  1976v  '  [based  on  the  original  design  by  Glick  In  1955] 


After  some  modifications  by  Vaughn  in  his  soot  studies  from  benzene 

other  authors  have  performed  studies  using  this  shock  tube  on  benzene 

( 36  37  3ft) 
and  cyclohexane  for  soot  formation.    '   '     For  more  information 

about  shock  tube  design  and  theory,  the  reader  should  consult  Gaydon  and 

Hurle.(34) 

The  shock  tube  has  been  built  using  type  304  stainless  steel  with 
inside  diameter  5.08  cm  and  wall  thickness  of  0.635  cm  capable  of 
withstanding  a  pressure  of  205  atm.   It  is  divided  into  two  pressure 
regions:  a  high  pressure  region  or  driver  section,  and  a  low  pressure 
region  or  experimental  section  where  the  test  section  is  located  (see 
Fig.  4). 

The  driver  section  is  of  variable  length  (2-3  m)  to  change  the 
reaction  time  in  the  test  section.   This  section  has  two  pressure  gauges 
to  measure  the  initial  driver  pressure  P.,  and  the  final  pressure  P, 
after  the  shock  wave  is  trapped  in  the  dump  tank  (y   3  atm) .   Helium  was 
selected  for  driver  gas  because  it  is  an  inert  gas  and  its  low  molecular 
weight  permits  the  strongest  shocks. 

The  experimental  section  includes  the  dump  tank  and  the  test 
section.   In  the  test  section  the  reactant  mixture  is  injected  by  means 
of  a  gas  line  connected  at  the  end  wall  (see  Fig.  5) .   Two  pressure 
transducers  are  located  in  this  section  19.75  cm  apart  in  order  to 
record  the  time  the  shock  wave  travels  this  distance  for  computing  other 
shock  wave  parameters  (see  discussion  in  next  section) .   A  ball  valve 
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Fig.  3.   Behavior  of  Shock  Wave  Parameters  (Time,  Pressure,  and 
Temperature)  After  Bursting  of  Diaphragm. 

a)  Conventional  Shock  Tube 

b)  An  (x-t)  diagram  showing  progress  of  the  shock  wave, 
the  rarefaction  fan  and  the  contact  surface  separating 
driver  and  experimental  gases. 

c)  The  pressure  distribution  along  the  tube  at  time  t. . 

d)  The  temperature  distribution  at  time  t  . 
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Fig.    4.      Shock  Tube  Diagram. 
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separates  this  region  from  the  rest  of  the  shock  tube.   At  the  other  end 
of  the  experimental  section  and  near  the  diaphragm  assembly  there  is  a 
40  liter  dump  tank  connected  in  an  obtuse  angle  to  the  shock  tube  which 
serves  as  the  shock  wave  trap.   The  experimental  or  driven  section  is 
5.9  m  long  and  argon  is  used  because  it  is  an  inert  gas  and  has  a 
constant  specific  heat  ratio. 

The  mechanism  for  bursting  the  aluminum  diaphragms  (OCLAD  0.04  cm) 
was  modified.   Originally,  the  diaphragms  were  ruptured  manually  with  a 
plunger  located  in  an  acute  angle  with  the  driver  section  near  the 
diaphragm  assembly.   A  new  system  described  in  detail  in  Appendix  B  was 
designed  by  the  author  in  conjunction  with  Kent  Sutton  (conducting 
experiments  on  benzene  pyrolysis) .   This  new  system  allows  a  better 
reproduction  of  the  experimental  conditions.   The  system  is  electrically 
driven  and  performs  by  de-energizing  an  electro-magnet  that  holds  the 
plunger  in  position  by  compressing  a  stainless  steel  spring.   Once  the 
magnet  is  turned  off,  the  spring  expands  and  drives  the  plunger  into  the 
diaphragm  assembly  and  the  shock  wave  is  initiated. 

Long  before  the  experiment,  a  mixture  bottle  (stainless  steel,  1 
liter)  is  prepared  with  calculated  amounts  of  toluene  and  argon  (see 
Appendix  A) .   The  concentration  of  toluene  is  then  checked  by  gas 
chromatography.   This  mixture  is  introduced  into  the  previously 
evacuated  test  section.   Pure  argon  at  the  same  pressure  is  then 
injected  into  the  rest  of  the  experimental  section  (the  test  section  is 
separated  from  the  rest  of  the  shock  tube  by  closing  ball  valve) .   Once 
these  two  steps  are  performed  the  gauges  (speedivac  0-760  torr)  that 
measure  the  pressure  in  the  experimental  section  are  shut.   At  the  same 
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Fig.  5.   Close  up  of  Test  Section  Showing  the  System  for 
Gas  Sampling. 
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time  that  the  first  steps  are  performed,  the  driver  section  is  also 
evacuated.   Then,  the  helium  driver  gas  is  introduced  to  a  prescribed 
pressure  and  the  shock  wave  is  ready  to  be  fired.   The  steps  taken  to 
start  the  shock  wave  are  performed  in  a  very  short  time;  these  are: 
swinging  the  ball  valve  in  line,  turning  off  the  electro-magnet,  and 
shutting  the  ball  valve  closed  again.   After  a  period  of  three  hours, 
the  samples  are  collected  from  the  test  section  and  analyzed  by 
chromatography.   A  more  detailed  explanation  of  the  shock-wave 
initiation  procedure  is  presented  in  Appendix  C. 

2.4.3  Determination  of  Shock  Wave  Parameters 

(341 
According  to  Gaydon  and  Hurle     when  a  compression  wave  is  formed 

in  the  low  pressure  zone  the  speed  of  sound  in  argon  is  given  by 

a  =  /yRTTTMW  (2.2) 

where  y  Is  the  specific  heat  ratio,  R  is  the  universal  gas  constant,  T 
is  the  temperature  before  the  shock  wave  is  fired  and  MW  is  the 
molecular  weight  of  argon.   For  a  temperature  of  T  =  298  K,  a.  is  321 
m/sec.   Then,  the  application  of  the  conservation  of  mass,  momentum,  and 
energy  equations  when  ideal  behavior  and  constant  specific  heat  are 
assumed,  yields  the  following  relations  used  to  calculate  the 
temperature  and  pressure  at  the  reflected  region  5  as 

T  .  T   [2(T-1)M2  +  (3-y)][(3Y-1)M2  -  2(Y-D]  ,,  ,. 

5     1  2   2  *■*■•->) 

5    L  (Y+l)  m2 


2YM2  -  (y-1)    (3Y-1)M2  -  2(y-1)| 

5  '1  Y"    11  r-2*2  J 


(2.4) 


and  M  -  u/a  (2.4.1) 

u  =  d/t  (2.4.2) 

where  M  Is  the  mach  number  of  the  incident  shock  wave ,  y   is  the  specific 
heat  ratio  of  the  gas,  d  is  the  distance  between  pressure  transducers 
(19.75  cm),  t  is  the  time  it  takes  for  the  shock  wave  to  travel  distance 
d,  and  T. ,  P.  are  the  temperature  and  pressure  in  region  1.   The  mach 
number  is  typically  2,  and  can  be  calculated  from  the  recorded  time  't' 
in  the  Fluke  (model  1952B)  time  counter.   The  shock  wave  passes  by  the 
two  pressure  transducers  installed  19.75  cm  apart  in  the  test  section, 
and  the  elapsed  time  is  recorded.   This  amplified  voltage  goes  also  into 
the  Tektronix  oscilloscope  where  a  pressure  history  is  recorded  which 
will  be  useful  in  calculating  the  cooling  rate  addition  to  the  reaction 
time  (see  Appendix  D  for  detailed  discussion) . 

2.4.4  Temperature  Verification  from  Reflected  Shock  Study 

Studies  by  previous  authors  using  the  KSU  shock 

.  .  (2,35,36,37,38)  .     .  .    .   ,    .        ,  „.    .   , 

tube  have  determined  an  increase  in  the  shock  wave 

temperature  and  pressure  when  reflection  is  allowed  to  occur  at  the  end 

wall.   A  study  for  comparing  the  ideal  temperature  and  pressure  (Eq. 

(2.3)  and  (2.4))  with  the  actual  values  in  the  shock  tube  was  performed 

due  to  the  importance  of  these  two  parameters  in  the  experiments. 

A  series  of  shock  experiments  were  performed  where  the  incident  and 

reflected  velocities  were  measured.   The  time  't  '  travelled  for  the 

incident  shock  was  measured  with  the  Fluke  timer  and  the  time  't  '  for 

the  reflected  shock  was  measured  by  using  the  oscilloscope.   These  two 

measured  parameters  't  '  and  't  '  were  necessary  to  compute  the  incident 
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and  reflected  mach  numbers  (M  ,  MO  by  letting  t  or  t„  =  t  In  Eq. 
(2.4.2)  and  the  Mach  number  M^  or  M  =  M  in  Eq.  (2.4.1)  where  d  was  kept 
constant  at  19.75  cm.  and  these  were  used  in  the  equation 


,     lV-(¥W 


|Y-r   Mi 


«H 


(2.5) 


where  M  =  mach  number  MT  or  Mp,  and  also  in  Eq.  (2.3)  for  T,.. 

For  a  temperature  T  =  298  K,  Eq.  (2.5)  was  used  to  compute  T  by 
letting  M  =  K   ,   T   -  T_,  and  T  =  T,.      Then,  the  ratios  between  T  /T 
and  T'/T.  were  set  equal  at  the  end  wall  as 

T    T  ' 
2    5 

1    2 


using  Eq.  (2.5)  again  for  the  ratio  T'/T  where  M  =  M_(tR) ,  the 
reflected  temperature  T'  was  obtained.   Collected  data  for  this  study  is 
shown  in  Table  4  for  long  and  short  shock  tube  geometries  using  argon  as 
the  driven  gas. 

The  calculated  Incident  and  reflected  shock  temperatures  (T,,  T  ') 
were  fitted  using  least  squares  methods  to  yield  the  following  empirical 
relationships: 


T5'  =  (0.9134  T5  +  15.8320)K   r2  =  0.9876  (2.7) 

800  <  T  <  1250  K 
long  geometry. 

and 
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Table  4.   Experimental  Data  for 

Temperature  Calibration 
of  the  Short  and  Long 
Shock  Tube  Geometries 


R 
(Usee) 


I 
(Usee) 


5 
(K) 


5 

(K) 


Short  Tube  Data 


500 

264 

1266 

1348 

502 

262 

1278 

1367 

517 

274 

1175 

1257 

499 

264 

1267 

1348 

536 

282 

1105 

1192 

529 

280 

1125 

1208 

497 

260 

1294 

1380 

493 

256 

1331 

1421 

532 

284 

1097 

1176 

499 

271 

1220 

1286 

294 

265 

1264 

1335 

507 

274 

1188 

1257 

479 

264 

1289 

1341 

492 

266 

1262 

1329 

493 

269 

1238 

1298 

Long  Tube  Data 


542 

273 

1157 

1269 

528 

370 

1192 

1295 

547 

289 

1058 

1145 

512 

275 

1178 

1252 

533 

279 

1129 

1220 

523 

280 

1135 

1212 

543 

287 

1073 

1160 

589 

316 

891 

976 

551 

291 

1043 

1131 

564 

312 

929 

999 

552 

310 

949 

1010 

541 

298 

1017 

1084 
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T5'  =  (0.9753  T5  -  45.4235)K   r2  =  0.9829  (2.8) 

1250  <  T  <  1700  K 
short  geometry. 


Equation  (2.3)  overpredicts  the  temperature  value  by  only  1%    (70  K) 
in  the  temperature  range  800-1250  K  and  by  5.5%  (equivalent  to  83  K)  in 
the  temperature  range  1250-1700  K.   Equations  (2.7)  and  (2.8)  were  used 
in  calculating  the  actual  shock  tube  temperature  from  experimental  data 
used  in  Eq.  (2.3)  to  calculate  T,. 

2.5  Post-Shock  Sampling  Systems 

Once  a  shock  wave  was  performed,  the  solid  and  gaseous  samples  were 
allowed  to  equilibrate  for  three  hours.   After  this  step,  the  samples 
were  collected  in  two  different  ways  depending  on  their  nature  (solid  or 
gaseous) . 

2.5.1  Gaseous  Collection 

The  gaseous  collection  is  the  simplest  way  of  sampling,  and  it  was 
performed  by  filling  two  75  cm3  stainless  steel  bottles  which  were 
previously  evacuated  to  1.32  x  10   atm. ,  provided  with  a  quick 
connection  in  one  side  and  an  injection  port  in  the  other  side.   The  two 
bottles  were  connected  to  the  end  wall  of  the  shock  tube  by  means  of  a 
three-way  valve  (see  Fig.  5).   The  first  bottle  (dump)  was  neglected  for 
analysis  and  was  used  only  as  a  storage  for  all  impurities  that  were 
present  in  the  sampling  lines.   These  impurities  came  in  the  form  of 
toluene  and  argon  from  the  mixture  bottle  which  was  originally  connected 
to  the  line  where  the  dump  bottle  was  also  connected.   Also,  air  could 
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have  leaked  into  the  lines  when  the  mixture  bottle  was  disconnected  from 
the  test  section. 

2.5.2  Solid  Collection 

The  high  molecular  products  (those  with  molecular  weight  greater 
than  that  of  toluene)  were  collected  in  a  cylindrical  and  removable 
aluminum  liner  located  in  the  test  section  of  the  shock  tube.   After 
each  experiment,  the  liner  was  immersed  into  a  solvent  (CS.)  for  at 
least  8  hours.   Then,  the  solvent  containing  the  solid  samples  was 
distilled  at  319  K;  the  volume  of  solution  was  reduced  to  1  cm3  and  kept 
for  further  analysis  in  the  Tracor  GC  analyzer.   The  liner  was  heated 
for  12  hours  to  evaporate  residues  of  CS.  and  other  impurities,  before 
it  was  used  again  in  a  new  experiment. 

Standard  samples  for  the  Tracor  chromatograph  were  prepared  using 

bibenzyl,  diphenyl  methane  and  biphenyl  in  solution  with  carbon 

—7    —8       —9 
disulfide  (10   ,  10   and  10   mol/cm3  respectively) .   Other  compounds 

present  in  the  samples  were  detected  by  comparing  retention  times  with 

known  standards. 


39 


3.0  RESULTS  AND  DISCUSSION 
The  pyrolysis  of  vapor  phase  toluene  was  studied  in  the  temperature 
range  800  -  1700  K  for  two  concentrations.   A  reaction  time  study  was 
performed  at  1283  K  in  order  to  measure  the  extent  of  toluene 
decomposition  as  a  function  of  reaction  time.   The  results  from  these 
experiments  can  be  divided  into  three  categories:   1)  reaction  time 
study,  2)  effect  of  temperature  on  toluene  decomposition,  and  3) 
modeling  of  toluene  decomposition. 

3.1  Reaction  Time  Study 

A  series  of  shock  tube  experiments  were  performed  where  the  length 
of  the  driver  section  in  the  shock  tube  was  changed  from  2  m  to  3  m  in 
order  to  change  the  reaction  time. 

As  it  was  discussed  in  the  previous  chapter,  when  the  diaphragm  is 
burst  a  rarefaction  fan  is  formed  in  the  driver  section.   This 
rarefaction  fan  is  reflected  from  the  end  wall  of  the  driver  section  and 
when  it  arrives  at  the  experimental  section,  the  reaction  is  slowed  down 
by  decreasing  the  pressure  and  temperature  in  the  test  section.   By 
changing  the  length  of  the  driver  section,  the  time  of  arrival  of  the 
cooling  fan  is  also  changed,  and  hence,  the  reaction  time  will  also  be 
changed.   The  oscilloscope  trace  (see  Fig.  6)  provides  the  shock 
pressure  history  where  the  flat  portion  represents  the  period  of  time 
when  pressure  and  temperature  P  and  T  are  constant  (uncorrected 
reaction  time).   Then,  the  inclined  section  is  where  the  pressure  and 

the  temperature  decrease  in  an  adiabatic  process.   The  temperature  can 

(2) 
be  calculated  from  the  following  equation, 
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Fig.  6.   Pressure  History  of  the  Shock  Wave  Recorded  by 
the  Tektronix  Oscilloscope. 

y-axis  is  pressure  transducer  response  in  volts. 

x-axis  is  time. 
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where 


(y-D/y 

T(t)  =  T5  (-£Mj        ,  (3.1) 


P(t)  -  Pj  -  [dP/dt]t  .  (3.2) 


The  time  derivative  for  the  pressure  dP/dt  can  be  measured  from  the 
oscilloscope  trace,  y   is  the  specific  heat  ratio  which  for  argon  and 
helium  is  5/3,  't'  is  the  time  it  takes  for  the  temperature  to  go  from 
its  highest  value  (T  )  to  room  temperature. 

The  conditions  in  these  experiments  were  as  follows 

Temperature:    1283  ±  35  K 
Average  Toluene  mol  fraction  in  Argon:   6.7  x  10~ 
Reaction  time:    0.8  -  3.0  msec 

The  concentration  of  the  unreacted  toluene  was  measured.   The 
general  form  for  the  integrated  expression  for  toluene  disappearance  can 
be  written  either  as 


te  ■[!]--  -Kxt     n  =  1  (3.3) 


«  [T]1""  -  [T]^"n  =  -  Knt(l-n)     n  4   1       (3.4) 

where  K  is  the  reaction  rate  constant,  [T]  is  the  concentration  of 

toluene  at  time  't'  and  [T]   is  the  initial  concentration  of  toluene, 
o 

Calculations  were  made  of  each  condition  by  use  of  Eqs.  (3.3)  and  (3.4) 

with  assumed  reaction  orders  of  1,  1.5,  and  2.   The  results  are  shown  in 

2 
Table  5  where  r  is  the  correlation  coefficient  for  the  least  squares 
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fit.   Least-squares  regression  analysis  of  the  data  resulting  from 
plotting  the  left  hand  side  of  the  equations  (3.3)  and  (3.4)  versus  time 
't'  for  each  assumed  order  are  shown  in  Tables  5  and  6,  and  it  indicates 
that  the  best  fit  to  the  data  is  obtained  when  the  reaction  order  of  one 
is  assumed.   The  correlation  coefficient  r  draws  this  conclusion  where 
a  value  close  to  one  indicates  a  better  fit.   The  first  order  plot  is 


shown  in  Fig.  7.   These  results  are  in  good  agreement  with  results 

initially  reported  by  Swarc 

authors.(13>17'19'20>22'31> 


(39) 
initially  reported  by  Swarc     and  later  by  other 


3.2  Effect  of  Temperature  on  Toluene  Decomposition 

A  series  of  experiments  were  performed  where  the  product  yields 

were  determined  as  a  function  of  temperature  for  two  initial 

—3  ? 

mol  fractions  of  toluene  in  argon:  7.4  x  10   and  1.25  x  10   .   These 

concentrations  were  selected  over  higher  values  so  that  a  constant 

specific  heat  ratio  (y)  for  the  mixture  could  be  assumed.   A  higher 

concentration  would  cause  an  undesirable  specific  heat  ratio  temperature 

dependence.   The  reaction  (corrected)  time  was  kept  constant  at 

1.97  ±  0.13  msec  for  the  temperature  range  1100  -  1700  K,  and  also  at 

3.25  ±  0.24  msec  for  the  range  800  -  1100  K. 

A  list  of  the  observed  products  in  these  experiments  is  shown  in 
Table  7.   These  compounds  were  detected  by  the  chromatographic  analysis 
comparing  retention  times  with  known  samples,  and  their  concentrations 
were  determined  by  comparison  with  calibration  standards. 

Plots  of  the  major  percent  product  yields  per  millisecond  of 
reaction  versus  reaction  temperature  is  shown  in  Figs.  8(a)  through  (h) . 
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Table  5.  Least  Squares  Fit  Results  for  Reaction 
Time  Study. 


Reaction  Order       K  (mole, K, sec) 

1.0 
1.5 

2.0 


2 
r  =  correlation  coefficient  for  Least  Squares  Fit. 

K  ■  Reaction  Rate  in  mole, cm,  sec  units. 


418.1 

0.8532 

1.45  x  106 

0.8162 

5.01  x  109 

0.7757 
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Fig.  7.   Toluene  Thermal  Decomposition  Curve  That  Verifies  First  Order 
Behavior.   T  =  1283  K,  Initial  Toluene  Concentration  ■  0.67% 


f°U3/ci:]  vl-q 
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Table  6.   Experimental  Data  used  for  Least  Squares  Analysis  in 
Reaction  Time  Study. 


t  (msec)    -£n[[T]/[T]  ]    2(l/[T]!s  -  1/IT]*5)     (1/[T]  -  1/[T]  ) 
o  o  o 


1.74 

0.456 

1.40 

* 

103 

4.28 

*  106 

1.55 

0.304 

8.74 

* 

102 

2.53 

*  106 

2.28 

0.877 

3.07 

* 

103 

1.08 

*  107 

0.92 

0.285 

9.17 

* 

102 

2.99 

*  106 

1.30 

0.301 

9.77 

* 

102 

3.19 

*  106 

1.23 

0.139 

4.27 

* 

102 

1.39 

*  106 

0.97 

0.30 

8.70 

* 

101 

2.51 

*  105 

48 


Table  7.   Products  from  Toluene  Decomposition 
Observed  in  this  Study. 


Methane 

CH4 

Ethylene 

C2H4 

Acetylene 

C2H2 

Propadiene 

C3H4 

Propyne 

C3H4 

Propene 

C3H6 

Vinylacetylene 

C4H4 

Diacetylene 

C4H2 

1,3  Butadiene 

C4H6 

Styrene 

C8H8 

Benzene 

C6H6 

Biphenyl 

C12H10 

Diphenyl  Methane 

C13H12 

Bibenzyl 

C14HH 

Also  traces  of  C„,  C. 

.„,  C.„  unidentified  compounds. 
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These  major  product  yields  are  acetylene,  benzene,  methane,  and 
ethylene,  and  their  percent  yield  is  calculated  by  dividing  the 
composition  of  the  product  by  the  initial  composition  of  toluene 
injected  into  the  test  section.   The  yields  of  minor  products  (C,,  C,, 
and  high  molecular  weight  hydrocarbons)  are  presented  in  Appendix  H. 
The  major  conclusion  from  these  plots  is  that  benzene  is  the  major 
product  formed  at  low  temperatures  reaching  a  maximum  yield  at  1425  K. 
This  is  the  first  time  that  benzene  has  been  observed  with  this  high 
yield  in  shock  tube  low  temperature,  thermal  decomposition  of  toluene. 
R.  Smith   '   '    observed  benzene  and  other  hydrocarbons  using  a 
Knudsen  cell  and  the  following  experimental  conditions: 


Temperature  Range:   1173  -  2073  K 
Reaction  Pressures:  9.87  x  10~  -  g.87  x  10   atm. 
Reaction  Time:   3  x  10   (M/TKsec 

M  =  molecular  weight 
T  -  temperature 
Range:  0.840  -  0.632  msec 


In  this  study,  higher  reaction  pressures  were  used  in  excess  of  10  atm 
the  reaction  time  was  maintained  at  around  2  msec  for  shock  waves  with 
temperatures  above  1100  K  and  at  around  3.25  msec  at  lower  temperatures. 
It  is  important  to  note  that  the  low  weight  product  yields  in  Smith's 
work  are  very  similar  to  the  ones  found  in  this  work  (see  Table  8)  with 
the  only  difference  in  the  yield  vs.  temperature  distribution,  where  the 
maximum  yield  of  most  species  occurs  at  1673  K  which  is  248  K  higher 
than  the  temperature  where  most  of  the  species  observed  in  this  study 
reach  their  highest  yield  (1425  K) .   However,  acetylene  is  first 
observed  at  1273  K  by  Smith  which  is  close  to  the  threshold  temperature 
for  the  production  of  such  species  in  this  study  (1300  K) . 
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Fig.  8.   Product  Yield  vs.  Temperature  Distribution  for  Major 
Products  in  Toluene  Pyrolysis. 

(a)   Acetylene,  Initial  Toluene  Mol  Concentration 
(ITC)  -  1.25% 


(b) 

Acetylen 

,  ITMF 

=  0.74% 

(c) 

Benzene, 

ITMF  = 

1.25% 

(d) 

Benzene, 

ITMF  = 

0.74% 

(e) 

Methane, 

ITMF  = 

1.25% 

(f) 

Methane, 

ITMF  = 

0.74% 

(g) 

Ethylene 

,  ITMF  i 

=  1.25% 

00 

Ethylene 

,  ITMF  ' 

=  0.74% 

(1) 

Toluene, 

ITMF  = 

1.25% 

(J) 

Toluene, 

ITMF  = 

0.74% 

Note:   solid  lines  represent  computer  model  results. 
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Table 

8.   Product  Yield1  ' 

-Temper 

ature  Distribution 

Results  f 

rom  R. 

Smith  Studies 

Temperature  (K) 

Species 

1173 

1273 

1373 

1473 

1573 

1673 

1773 

1873 

2073 

CH 

CH4 
C2H2 
C2H3 
C2H4 
C3H3 

0.31 

1.5 

11.0 

24.0 

15.0 

5.2 

2.1 

0.26 

0.44 

2.0 

4.9 

4.9 

5.0 

3.8 

0.4 

0.02 

0.09 

2.0 

16.0 

37.0 

65.0 

78.0 

80.0 

0.14 

0.78 

1.7 

3.0 

4.0 

4.5 

0.12 

0.09 

0.20 

5.2 

3.2 

2.5 

1.5 

i-O.Ol 

0.09 

1.3 

5.0 

5.4 

1.2 

0.4 

0.2 

C3H4 

C3H5 

C4H2 

C4H3 
C?Hf 

CH 

C5H4 
C5H5 
SH6 

0.19 

0.78 

0.8 

0.3 

0.3 

0.5 

0.06 

0.21 

0.25 

0.1 

•vO.01 

0.01 

0.68 

4.1 

6.3 

5.5 

5.0 

12.0 

0.03 

0.20 

0.4 

0.7 

0.4 

0.3 

0. 

0.03 

0.20 

0.8 

1.2 

0.7 

0.3 

0.02 

0.33 

0.8 

0.7 

0.1 

0.01 

1.10 

0.3 

0.2 

0.01 

0.18 

1.2 

3.5 

2.8 

1.0 

0.1 

0.08 

0.13 

0.31 

0.3 

0.1 

C6H2 

0.01 

0.2 

0.6 

0.35 

0.3 

0.3 

C6H5 

C6H6 
CH 

0.24 

0.7 

0.8 

0.36 

0.1 

0.06 

0.15 

0.64 

1.1 

1.3 

0.7 

0.25 

0.13 

0.90 

4.4 

10.0 

11.0 

7.0 

3.0 

0.5 

0.03 

0.17 

1.2 

1.7 

0.6 

0.05 

C7H6 

1.0 

1.0 

0.3 

0.01 

CH7 

14.0 

30.0 

40.0 

37.0 

16.0 

2.5 

0.1 

0.02 

86.0 

69.0 

56.0 

33.0 

8.4 

1.1 

C8H2 

0.03 

0.01 

C8H6 

0.07 

0.45 

0.6 

0.2 

0.4 

Vs 

0.10 

0.33 

0.15 

C9H7 

0.08 

0.5 

0.4 

0.06 

ClX 

0.01 

0.07 

0.08 

0.03 

C10H8 

0.02 

0.2 

0.4 

0.15 

0.01 

C  1H7 

0.01 

0.09 

0.06 

C2H8 

0.03 

0.18 

0.2 

0.05 

0.01 

C12H9 

0.20 
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0.01 
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The  yield  of  methane  as  a  function  of  temperature  is  shown  in  Figs. 
8(e)  and  (f ) ,  and  it  reaches  a  maximum  value  at  1500  K.   Methane  first 
appears  at  1300  K  which  is  somewhat  higher  than  the  temperature  at  which 
benzene  begins  to  appear.   Smith  observed  methane  reaching  its  highest 
yield  at  1773  K  and  appearing  first  at  1373  K.  Ethylene  first  appears 
at  a  lower  temperature  (1250  K)  and  reaches  a  maximum  at  1425  K.   The 
yields  of  ethylene  as  a  function  of  temperature  are  shown  in  Fig.  8(g) 
and  (h). 

Acetylene  is  observed  at  1300  K  and  increases  rapidly  above  1425K 
when  the  yield  of  benzene  and  ethylene  start  decreasing.   This  is  shown 
in  Figs.  8(a)  and  (b) .  This  rapid  increase  in  the  yield  of  acetylene  as 
a  function  of  temperature  coupled  with  the  rapid  decrease  in  toluene 
concentration  and  the  yield  of  benzene  suggests  that  fragmentation  of 
the  tolyl  radical  becomes  more  important  with  increasing  temperature. 

The  yields  of  the  high  molecular  weight  products,  biphenyl 
(C.  H  ),  diphenylmethane  (C..H  .)  and  blbenzyl  (C..H.,)  are  orders  of 
magnitude  smaller  than  the  yields  of  benzene  and  these  are  shown  in 
Appendix  H.   These  results  indicate  that  the  cleavage  of  the  methyl 
group  from  toluene, 

C6H5-CH3  — *  CgHy  +  CH3'  (3.5) 

does  not  compete  with  the  C-H  bond  cleavage 

C6H5-CH3  — »  C6H5CH2'  +  H-  .  (3.6) 

A  series  of  experiments  were  performed  in  which  the  toluene 
mol  fraction  in  the  bath  was  maintained  constant  at  about  0.74  x  0.002% 
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and  the  corrected  reaction  time  was  kept  at  2.01  ±  0.21  msec.   The 
temperature  was  varied  from  1100  to  1700  K.   The  resultant  major  product 
yield  curves  are  shown  in  Figs.  8(b) , (d) , (f ) , (h)  and  the  other  product 
yield  curves  can  be  found  in  Appendix  H.   By  comparing  the  curves  from 
the  previous  experiment  with  the  new  ones,  a  similar  product-  temperature 
profile  is  observed  which  means  that  the  decomposition  of  toluene  is  of 

first  order.   This  is  because  when  the  original  initial  toluene  mol 

-2  -3 

fraction  of  1.25  x  10   was  decreased  by  half,  (to  around  7.4  x  10  ), 

the  product  yields  also  decreased  by  half.   Also,  the  toluene 

-2 
decomposition  curve  is  very  similar  to  the  one  from  1.25  x  10   toluene 

mol  fraction  in  argon  (see  Fig.  8(i),(j)). 

3.3  Modeling  of  Toluene  Decomposition 

The  present  section  will  discuss  the  reaction  mechanisms  involved 
in  the  thermal  decomposition  of  toluene.   First,  a  small  review  of  the 
literature  describing  some  of  the  proposed  mechanisms  in  toluene 
pyrolysis.  Then,  an  explanation  of  the  modeling  techniques  used  to 
describe  the  experimental  data  from  the  present  work. 

3.3.1  Proposed  Reaction  Mechanisms 

A  review  of  the  literature  shows  an  agreement  among  the  toluene 
pyrolysis  experimentalists  in  the  first  steps  involved  in  toluene 
decomposition.   No  matter  how  the  pyrolysis  is  performed  (flow  system, 
static  system  or  shock  tube)  or  what  the  conditions  of  concentration, 
pressure  and  temperature  are,  the  first  three  proposed  reactions  are 

CyH8  ►  C7Ry  +  H-  (3.5) 
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H-  +  C7Hg  ►  C?H7-  +  H2-  (3.6) 

H*  +  C7H8 C6H6  +  CH3-  °-7) 

(39) 
In  1947,  M.  Swarc     was  the  first  experimentalist  to  propose  a 

reaction  mechanism.   Working  with  a  flow  reactor  he  observed  hydrogen, 

methane  and  bibenzyl  as  the  relevant  products.   Based  on  these 

observations  he  proposed  a  mechanism  that  included  Eqs.  (3.5),  (3.6), 

(3.7)  and  also  the  following  for  bibenzyl  formation 

CH3-  +  C?H8  ►  C?H7'  +  CH4  ,  (3.8) 

2  CyH7-  ►  C14H14  (bibenzyl)  .      (3.9) 

He  also  reported  a  reaction  rate  constant  for  reaction  (3.5)  given  by 

Kj  =  2.1  x  1013  exp(-324.5  (kJ) /RT)  sec"1 

Following  the  work  performed  by  Swarc,  Steacie,  et  al.     in  1954 
argued  about  the  possibility  of  the  reaction 

C7H8  ►  C6H5-  f  CH3-  (3.10) 

to  explain  the  formation  of  methane,  based  on  the  fact  that  they  found  a 
higher  activation  energy  for  reaction  (3.5)  (376.8  kj)  than  the  value 
reported  by  Swarc. 

Later,  M.  Takahasi(18'41'  in  1959  and  S.  J.  Price(19)  in  1961 
working  with  flow  systems,  proposed  a  new  set  of  reactions.   Takahasi 
explained  the  formation  of  hydrogen  and  methane  by  hydrogen  atom  and 
methyl  radical  abstraction 


^  C?H8  +  H "  C7H7-  +  H2  '  "•« 
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CyH8  +  CH3-  ►  C7Hy  +  CH4  .  (3.11) 

Price  also  included  reactions  (3.6)  and  (3.11)  into  the  mechanism  and 
added  the  reaction 


2  CH3-  ►  C2H6  ,  (3.12) 

to  explain  the  formation  of  ethane 

In  1971,  Peacock  et  al.,     introduced  two  new  reactions  to  the 
mechanism 


CH3-  +  H2  ►  CH4  +  H-  (3.13) 

C7H?-  +  H2  ►  C?H8  +  H«  (3.14) 

Reaction  (3.13)  produces  more  hydrogen  atoms  which  are  used  by  reaction 

(3.7)  to  form  benzene  and  reduces  the  formation  of  bibenzyl  via  reaction 

(3.8)  and  (3.9).  R.  Smith  using  Knudsen  Cells  (in  1979)  argued  that 
acetylene  was  the  main  precursor  for  soot  formation  based  on  his 
acetylene  yields  from  toluene  pyrolysis.  He  proposed  the  reactions: 

C7H?-  ►  C5H5-  +  C2H2  (3.15.a) 

C7H?-  ►  C4H4  ♦  C3H3-  (3.15.b) 

C5H5-  ►  C5H3-  +  H2  (3.16) 

C4H4  ►  C4H2  ♦  H2  (3.17) 

C4H4  "  C4H3  +  H> 

where  other  compounds  like  vinylacetylene  and  diacetylene  also  played  an 
important  role  in  soot  formation.  A  summary  of  the  initial  reactions  in 
toluene  pyrolysis  proposed  by  these  authors  is  shown  in  Table  9. 
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3.3.2  Modeling  of  Experimental  Data 

The  process  of  modeling  the  experimental  product  yields  is  a  very 
simple  one.   First,  all  the  possible  reactions  (Fig.  9)  and  other 
reactions  that  involve  benzene  decomposition  were  included  in  the 
pyrolysls  mechanism.  Reaction  rate  constant  values  were  provided  by  the 
literature.   A  system  of  first  order  differential  equations  was  set  up 
of  the  form 


£$■  =  K^AHCt    +  K^BJtA]    +  KjJCHB] 


^-   =  KjJB]  [A]  +  Kj  [C]  [A]  +  K±[B]  [A] 


where  [A],  [B],  [C],  etc,  are  the  concentration  of  each  species  present 
in  the  system,  K  values  are  the  rate  constant  for  reactions  i,  j,  k, 
etc.,  and  dt  is  the  time  derivative.   The  literature  values  for  the  rate 
constants  are  shown  in  Fig.  9. 

The  system  of  equations  was  solved  for  each  temperature  value  from 
1100  K  to  1700  K  by  using  a  computer  program  which  is  described  in 
Appendix  F.   This  computer  program  was  partially  provided  by  the  I. M.S. 
Library  routine  DGEAR  written  in  Fortran  (double  precision) .   The  method 
of  solution  for  this  system  of  equations  is  basically  the  backward 
differentiation  formula  also  called  Gear's  stiff  method.   The  computer 
programming  system  was  used  in  the  following  way.   First  of  all,  a 
system  of  first  order  differential  equations  was  set  up  including  the 
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Table  9.   Chemical  Reaction  Proposed 

by  the  Literature  in  Toluene 
Pyrolysis 


C7H8  *   C7H7-  +  H' 

H-  +  C?Hg  ►  C?H  •  +  H- 

H'  +  C7H8  h  C6H6  +  CH3' 

CH3-  +  C7Hg  ►  C7H?.  +  CH4 

2  C7H?.     ►  CUHU 

C7H8  "  C6H5  +  CH3- 

C7Hg  ♦  H ►  C7H;.  ♦  H2 

C7H8  +  CH3 ►  C7H?.  ♦  CH4 

2  CH3 ►  C2H6 

CH  •  +  H.  ►  CH.  +  H- 

3  2         4 

C7H7-  +  H2  ►  C?H8  ♦  H- 

2  c?H7  ►  C5H5  ♦  C2H2 

C5H5  "  C5V  +  H2 

2  C4H4  ►  C4H2  +  H2  ♦  C4H3  ♦  H- 
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Fig.  9.  Reported  Reactions  in  Toluene  Pyrolysis  and  their 
Respective  Reaction  Rate  Constants  used  Initially 
in  Computer  Modeling.   Units:  mole,  sec,  k,  cm3. 

(231 

References:   1)   Skinner  and  Subba  Rao 

2)  Empirical  Value         .... 

3)  Astholz,  Durant  and  Troe!:. 

4)  Kern,  Skinner,  et ,-al.    ' 

5)  Peacock,  et.  fl.»   ' 

6)  Asaba  &  Fujii1-  k) 


C?Hg  ►  C7Hy-+  H-  K  =  2.7xl04  exp(-347 .5(kJ) /RT) (1) 

H'+C  Hg  ►  C6Hfi+  CH3-  K  =  2.8xl013  exp(-15.5(kJ) /RT) (1) 

H-+C?H8  ►  C?H7-+  H2  K  =  7.6xlO~5T5-5  exp(-1424(J) /RT) (1' 

2  C7H7-  ►  C5H5'+  C2H2+  C4H4+  SV   K  =  8-3xl°l?  exp(-455(kJ) /RT) (3) 

H-+  C6H6  >•  CfiH  •  +  H2  K  =  101345  exp(-67(kJ)/RT)(4) 

CfiH6  >•  C4H4+  C2H2  K  =  1014-11  exp(-368(KF)/RT)(4) 

CfiH6  ►  C6H5-+H-  K  =  1015-08  exp(-343(kJ)/RT)(4) 

C6H5-  ►  C4H3-+  C2H2  K  =  1011-3  exp(-214(kJ)/RT)(4) 

C4H3-  1-  C4H2+  H-  R  -  2.82xl012  exp(-59(kJ)/RT)  (5) 

C?H7-+  H2  ►  C?H8+  H-  K  -  107-85  exp(-42(kJ)/RT)(6) 

C6H6  +  CfiH5 >•  C12H1Q+  H-  K  =  109-9  exp (-72 (kJ)  /RT)  (6) 

CH3-+  C6H6  ►  C7Hg+  H-  K  =  109"5  exp(-87(kJ)  /RT)  (6) 

CM4+  CfiH5-  ►  C?H8+  H-  K  =  107-2  exp(-39(kJ) /RT) (6) 

CH3-+  CgH6  >•  CH4+  C6H5.  K  =  lO11-1  exp(-224(kJ) /RT) (5) 

C?H8+  H2  ►  C6Hfi+  CH4  K  =  lO8-9  exp(-46(kJ)/RT) (6) 

CH4+  C5H5-  ►  CH3-+  C6Hg  K  =  109-6  exp(-34(kJ)/RT)  (6) 

CfiH5-+  H2  ►  C6Hg+  H«  K-  lO11-5  exp(-29(kJ)/RT)(6) 

C12H10+  H'  *  C6H6+  C6H5*  K  =  1()U'5  «P(-29(KJ)/RT)(6) 

C6H5-  ►  C2H2+  C4H2+  H-  K  =  1.2xl015  exp(-343(KJ) /RT) (1) 

CH3-+  C7Hg  ►  C?H7-+  CH4  K  =  lO10  exp(-15(KJ) /RT) (2) 
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stable  and  radical  species  shown  in  Fig.  10(1).   For  each  species,  the 
differential  equation  represented  the  rate  of  production  and 
dissappearance  of  such  species  observed  in  the  reaction  mechanism  first 
proposed  (Fig.  9).   Once  the  system  of  equations  was  set  up  for  each 
species,  literature  values  for  the  reaction  rate  constants  were  input  in 
the  program.   Also,  other  parameters  such  as  initial  toluene 
concentration,  reaction  time,  stepsize,  and  accuracy  were  input.   A 
value  of  2.87  x  10   mol/cm3  for  toluene  concentration  was  selected 
because  it  was  basically  the  value  used  in  the  experimental  data 
(equivalent  to  1.25  x  10   toluene  mol  fraction).   A  value  of  2  msec  was 
used  for  reaction  time  in  a  temperature  range  from  1100  K  to  1700  K; 
however,  the  product  yields  were  recorded  for  1  msec  of  reaction  time 
because  the  experimental  data  was  normalized  to  1  msec  except  for 
toluene  which  was  observed  at  2  msec.   Since  the  reaction  rate  constant 
values  were  temperature  dependent,  the  system  of  differential  equations 
had  to  be  solved  for  each  input  temperature  value. 

A  sensitivity  analysis  was  performed  with  the  initial  set  of 
equations  found  in  the  literature  (Fig.  10).   The  results  from  this 
analysis  generated  nine  reactions  (Fig.  11)  which  describe  the 
temperature  profiles  for  benzene,  methane,  and  acetylene  as  well  as  the 
toluene  decomposition  in  the  temperature  range  1100-1700  K.   Figures 
8(a) , (c) , (e)show  the  comparison  between  the  experimental  results  and  the 
expected  yields  from  the  model  where  the  solid  line  represents  the 
computer  values.   From  these  figures  it  is  observed  an  almost  perfect 
fit  for  benzene  and  acetylene  yields,  and  a  rather  close  fit  for  the 
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Fig.  10.   Species  used  for  Sensitivity  Analysis 

(1)  Species  used  initially 

(2)  Species  used  in  final  model 


Stable  Species 
(1)  (2) 


C7H8 

C7H8 

C6H6 

C6H6 

»2 

H2 

C2H2 

C2H2 

C4H4 

C4H4 

C4H2 

CH4 

C12H10 

CH. 
4 

Radicals 

C7H?. 

C7H7- 

H- 

H- 

CHj- 

CH3. 

C5H5- 

SV 

C3H3. 

C3H3- 

C6H5- 

C6H5- 

<w 

C4H3- 
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Fig.  11.   Final  Reactions  Used  for  Computer  Modeling  of 
Experimental  Data.   Units:  mole,  K,  sec,  cm3. 


K 
C?Hg  ►  C?H7+  H-  Kj  =  2.7xlOU  exp(-307(kJ)  /(RT) 

K2  13 

H-+  C?H8  ►  C6Hfi+  CH3-  K2  =  0.4x10  J  exp(-15(kJ) /RT) 

K 

H'+  C7H8  "■  C7H7'+  H2  K3  =  7.6xl0"5T5-5  exp(-1424(J) /RT) 

K 
2  C7H7"  *  C5H5*+  C2H2+  C4H4+  SV   K4  =  8-3xl°17  exp(-455(kJ)/RT) 

V 

H-+  C6H6   **  C6H5-+  H2  K5  =  1013'45  exp(-67(kJ)/RT) 

If 


C6H6  "   C4H4+  C2H2  h  '   10      exp(-318(kJ)/RT) 

K 
C6HS-  ►  C4H3-+  C2H2  K?  =  1015-08  exp(-343(kJ)/RT) 


C7H7'+  H2  "  C7H8+  H'  K8  =  2-82xl°12  exp(-59(kJ)/RT) 

CH3-f  C7Hg  — ^  C7H7-+  CH4  (a)K9  =  1013  exp(-67(kJ)/RT) 


(a) 

Value  generated  by  sensitivity  analysis. 
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methane  yields.   This  last  discrepancy  is  due  to  the  possible  existence 
of  other  reactions  which  were  not  included  in  the  reaction  mechanism, 
and  also  to  the  experimental  errors  associated  with  each  preexponential 
term  and  energy  of  activation  in  the  reaction  rate  constants  reported  in 
the  literature. 

It  is  important  to  note  that  the  originally  reported  values  for  the 
energy  of  activation  in  reaction  (1)  (Fig.  11)  (347.5  kj)  (19)  was 
changed  by  the  sensitivity  analysis  to  307,326  J  and  318,212  J.   It  is 

also  important  to  note  that  different  values  have  been  reported  for 

(  39) 
this  crucial  reaction.   For  example,  Swarc     reported  a  value  of 


324.5  kJ,  Takahasi^   '    a  value  of  308.6  kj  in  the  temperature 
1011-1122  K,  and  Steacie,  et  al.     reported  a  value  of  376.8  kJ. 


range 


Also,  the  preexponential  term  in  reaction  (2)  was  changed  from  0.4  x 

13  13 

10   to  2.8  x  10   ;  and  the  reaction  rate  constant  of  reaction  (9)  was 

generated  by  the  analysis.   This  last  value  is  an  empirical  result  used 

to  fit  the  experimental  data  and  should  not  be  interpreted  as  a 

calculated  or  measured  reaction  rate  value.   The  validity  of  these 

changes  is  supported  by  the  fact  that  the  literature  values  were 

obtained  in  different  experimental  conditions  to  those  in  the  present 

work. 

3.4  Summary  of  Results 

The  reaction  time  study  in  shock  tube  pyrolysis  of  vapor  phase 

toluene  led  to  the  conclusion  that  the  reaction  mechanism  is  of  first 

order.   This  conclusion  is  also  supported  by  the  concentration  studies 

where  the  toluene  decomposition  and  product  yield  results  at  two 
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_2 
different  initial  toluene  mol  fraction  in  argon  (1.25  x  10   and  7.4  x 

10  )  show  that  as  initial  concentration  is  increased  by  a  factor  of 

two,  the  product  yields  are  also  increased  the  same  factor. 

Product  yields  identified  by  gas  and  liquid  chromatography  were 

methane,  ethylene,  acetylene,  propadiene,  propyne,  propene, 

vinylacetylene,  diacetylene,  1-3  butadiene,  benzene,  and  high  molecular 

weight  compounds  such  as  bibenzyl  biphenyl,  diphenyl  methane,  styrene, 

and  hypothesized  traces  of  C  ,  C_,  C,„,  Cla  compounds. 

o    y  1U    lo 

The  observation  of  benzene  high  yields  at  low  temperature  shock 
tube  study  is  of  special  importance  because  it  is  the  first  time  that 
has  been  observed  in  such  conditions  and  initially  hypothesized  reaction 
mechanisms  for  bibenzyl,  biphenyl  and  diphenyl  methane  formation  seem  to 
play  a  smaller  role  in  toluene  pyro lysis.   Modeling  of  benzene,  methane 
and  acetylene  yields  show  a  very  good  agreement  between  experimental 
results  and  theory. 

The  overall  disappearance  rate  of  reaction  was  calculated  to  be 
K  -  2.5  x  10  exp(-107  (kJ) /RT)  sec"  for  the  temperature  range 
1100-1700  K. 
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4.0   CONCLUSIONS  AND  RECOMMENDATIONS  FOR  FURTHER  STUDY 

The  pyrolysis  of  vapor  phase  toluene  was  performed  in  the  Kansas 
State  University  Single  Pulse  Shock  Tube  for  three  sets  of  conditions: 
1)  reaction  time  study,  2)  effects  of  temperature  studies,  and  3) 
effects  of  concentration  studies. 

For  the  reaction  time  study,  where  the  temperature  and  initial 

toluene  mol  fraction  in  argon  were  kept  constant  at  1283  ±  35  K  and 

-3 
6,7  x  10   respectively,  and  the  reaction  time  was  varied  from  0.8  to 

3.0  msec,  the  order  of  the  reaction  for  toluene  decomposition  was 

determined  to  be  first. 

Also,  concentration  studies  showed  that  the  pyrolysis  of  toluene  is 

of  first  order  or  autocatalytic.   This  conclusion  has  been  drawn  after 

-3 
increasing  the  initial  toluene  mol  fraction  in  argon  from  7.4  x  10   to 

-2 
1.25  x  10   at  a  constant  reaction  time  of  2.01  ±  0.21  msec  in  the 

temperature  range  1100  -  1700  K.   The  results  showed  an  increase  by  a 

factor  of  2  of  all  the  species  yields  after  the  initial  toluene  mol 

-3  -2 

fraction  in  argon  was  doubled  from  7.4  x  10   to  1.25  x  10   ,  that  is 

explained  by  the  first  order  reaction  conclusion. 

The  main  product  yields  from  toluene  pyrolysis  were  acetylene, 

benzene,  methane,  and  ethylene.   Minor  products  were  propyne,  propene, 

propadiene,  diacetylene,  vinylacetylene,  styrene,  bibenzyl,  biphenyl  and 

diphenyl  methane.   The  yield  of  the  first  four  compounds  is  of  great 

importance  because  the  initial  hypothesis  that  motivated  the  present 

work  to  explain  the  formation  of  intermediate  species  to  yield  soot  was 

convincingly  rejected.   This  hypothesis  included  the  reactions 
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C7Hg  ►  C?H7  +  H-  (4.1) 

C7Hg  ►  CgH5.  +  CH3-  (4.2) 

leading  to  the  formation  of  bibenzyl,  biphenyl  and  diphenyl  methane. 
The  production  of  benzene  and  methane  in  large  amounts  suggests 
that  the  toluene  pyrolysis  at  these  conditions  follows  a  different 
reaction  pattern  which  is 

C7H8  ►  C?H7  +  H-  (4.3) 

H>  +  C7H8  "   C6H6  +  CH3'  <*•« 

CH3-  +  C7H8  "  SV  +  CH4  <-A-51 

Also,  the  increasing  production  of  acetylene  with  increasing 
temperature  which  starts  when  all  the  other  species  begin  to  decline  in 
production  yield,  suggests  that  the  fragmentation  process  becomes 
predominant  at  around  1400  K. 

The  use  of  computer  programs  for  toluene  decomposition  produced  a 
reaction  mechanism  which  gives  an  estimate  of  the  acetylene,  benzene  and 
methane  yields  as  well  as  the  toluene  decomposition  with  a  very  good 
accuracy. 

A  few  recommendations  for  further  studies  in  toluene  pyrolysis 
include  the  addition  of  some  more  sophisticated  analytical  devices  to 
measure  the  behavior  of  free  radicals  as  well  as  the  use  of  isotopically 
labelled  toluene  to  observe  the  production  pattern  of  certain  radicals 
which  are  important  in  the  proposal  of  reaction  mechanisms.   Also,  a 
qualitative  and  quantitative  analysis  should  be  performed  on  high 
molecular  products  to  observe  their  behavior  and  role  in  the  formation 
of  soot  at  higher  temperatures.   Finally,  mixture  studies  should  be 
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performed  where  other  compounds  could  be  added  to  toluene  and  then 
pyrolyzed  In  order  to  observe  the  suppression  of  acetylene  yield,  and 
hence  the  suppression  of  soot. 
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Appendix  A  -  Gaseous  Mixture  Preparation 

The  preparation  of  the  Toluene/Argon  mixture  to  be  used  in  the  test 
section  of  the  shock  tube  was  performed  in  the  following  way.   A  one 
liter  stainless  steel  bottle  was  equipped  with  an  injection  port  and  a 
quick  connector.   The  bottle  was  evacuated  and  liquid  toluene  was 
injected.  The  amount  of  toluene  and  argon  to  be  injected  was  determined 
from  the  ideal  gas  law 

PV  =  nRT  (A.l) 

where  P  is  pressure,  V  is  the  volume,  n  is  the  number  of  moles,  T  is 
temperature  and  R  is  the  universal  gas  constant.   In  Eq.  (A.l)  P  can  be 
written  as 

PB  =  P (toluene)  +  P (argon)  (A. 2) 

where  P (Toluene)  =  x  *  P„  (A. 3) 

X  is  the  mol  fraction  of  toluene  in  argon  in  the  mixture  bottle  which  in 

—2  T 

this  study  was  changed  from  1.25  x  10   to  7.4  x  10   .   The  vapor 

pressure  of  toluene  at  room  temperature  was  calculated  to  determine  the 

total  pressure  in  the  mixture  bottle.  The  temperature  dependence  of  the 

(3) 
toluene  vapor  pressure  is  reported  by  Hancock    as 

logtP(mmHG)]  -  6.95508  -  Wffifffaff|        (A.4) 

At  an  average  room  temperature  of  25°C  the  vapor  pressure  of  toluene  is 

approximately  28  mmHg.   Thus,  using  this  value  for  P(toluene)  in 

Eq.  (A. 3)  the  bottle  total  pressure  was  easy  to  determine  for  a 
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prescribed  initial  toluene  mol  fraction.   The  volume  amount  of  toluene 
(V  )  injected  was  then  determined  from  Eq.  (A.  1)  as 


v  _  P(toluene)  V(bottle)  MW(toluene) 
T  "        R  T  p (toluene)  (A,5) 


where  MW  is  the  molecular  weight. 

The  conditions  applied  in  the  present  study  were 


Toluene  mol  fraction  Toluene  Volume  Argon  Pressure 
(cm3) (atm) 

0.0125               0.163  3.69 

0.0074               0.098  3.04 


86 


Appendix  B  -  Design  of  the  New  Electromagnetically  Driven  Plunger 

The  KSU  Single  Pulse  Shock  Tube  was  provided  with  a  new  feature 
which  allowed  a  better  reproduction  of  the  shock  wave  temperature.   This 
new  feature  was  the  electromagnetlcally  driven  plunger  which  burst  the 
aluminum  diaphragm.   A  picture  of  the  new  device  is  shown  in  Fig.  B.l, 
where  two  stainless  steel  springs  can  be  seen  and  the  operating 
mechanism  is  similar  to  that  of  a  pinball  machine.   Spring  A  is  capable 
of  compressing  when  the  flat  end  of  the  plunger  is  held  by  the 
electromagnet;  when  this  magnet  is  de-energized,  the  plunger  is  released 
with  a  force  of  about  445  N  towards  the  diaphragm  assembly  in  the  shock 
tube,  rupturing  the  aluminum,  and  initiating  the  shock  wave.   Spring  B 
is  used  to  stop  the  plunger  so  that  it  produces  an  even  rupture  in  the 
diaphragm,  and  hence,  a  good  shock  structure. 

This  new  device  was  tested  for  a  number  of  shocks  run  at  a  pressure 
Pj^  of  0.47  atm  and  P  =■  19  atm  using  helium  and  argon  as  the  driver  and 
driven  gases  respectively.   Table  B.l  shows  the  experimental  data  for 
this  test,  where  the  average  temperature  was  1388  K  with  a  deviation  of 
±  37  K  or  2.3%.  This  result  is  much  better  than  shocks  run  with  the 
original  manually  driven  plunger  where  temperature  replications  had  a 
deviation  range  of  about  ±  10  to  ±  15%. 

The  deviation  observed  using  this  new  device  was  mostly  attributed 
to  the  presumed  difference  in  cold  working  the  material  when  it  was 
manufactured.   For  this  reason,  the  aluminum  diaphragms  were  annealed  at 
623  K  for  two  days  before  they  were  used.   It  was  also  established  that 
for  even  better  replication  of  shock  waves,  small  amounts  of  aluminum 
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diaphragms  should  be  annealed  because  when  they  are  used  2  or  3  days 
after  annealing  larger  deviations  in  temperature  occur. 


Fig.  B.l.  New  Electromagnetically  Driven  Plunger  in 
the  KSU  Shock  Tube. 

(A)  Driving  Spring 

(B)  Stopping  Spring 
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Table  1 

3.1.   Experimental  Data 

for  Tes 

ting  New 

Plunger 

Run 

At  (msec) 

T^  00 

1 

0.2609 

1378 

2 

0.2598 

1389 

3 

0.2603 

1384 

4 

0.2604 

1383 

5 

0.2558 

1429 

6 

0.2617 

1370 

7 

0.2549 

1438 

8 

0.2561 

1426 

9 

0.2639 

1350 

10 

0.2590 

1393 

11 

0.2612 

1372 

12 

0.2633 

1352 

91 


Appendix  C  -  Procedure  for  a  Shock  Tube  Experiment 

The  following  are  the  steps  to  perform  a  shock  tube  experiment 
after  a  sample  mixture  (in  this  case  toluene/argon  mixture)  is  prepared, 
and  the  aluminum  diaphragm  is  annealed  for  2  days  at  623  K.   The  steps 
can  be  followed  by  looking  at  Fig.  CI  where  a  schematic  picture  of  the 
shock  tube  is  shown. 

-  Insert  diaphragm  into  shock  tube  diaphragm  assembly  (Section  A) 
and  turn  on  plunger  magnet  (at  B) 

-  Evacuate  the  whole  shock  tube  to  1.32  x  10   atm 

-  Once  shock  tube  is  evacuated  close  valve  C,  so  the  experimental 
section  is  isolated.  Also  close  ball  valve  (D)  in  test  section. 

-  Fill  up  test  section  with  sample  mixture  to  desired  pressure  P.  . 

-  Fill  up  the  rest  of  the  experimental  section  with  pure  argon  to 
pressure  P  . 

-  Close  gauge  1  in  experimental  section. 

-  Close  valve  E  and  fill  up  driver  section  with  helium  to  desired 
pressure  P.. 

-  Once  shock  tube  is  ready,  close  all  valves  and  perform  the 
following  in  the  shortest  period  of  time  possible: 

1.  Open  ball  valve  D  in  the  test  section 

2.  Turn  off  the  magnet  (switch  B)  to  rupture  diaphragm 

3.  Close  ball  valve  D  again. 

-  Record  final  pressure  P5  in  gauge  4,  time  for  shock  wave  to  travel 

between  transducers  (at  F) ,  get  a  pressure  history  picture  from 

oscilloscope  (at  G) ,  and  then  wait  three  hours  before  collecting 

solid  and  gaseous  samples. 
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Fig.  C.l.   Schematic  Representation  of  the  Shock  Tube. 
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Appendix  D  -  Determination  of  Cooling  Rates 

The  determination  of  cooling  rates  in  the  present  study  was 

performed  following  the  procedure  described  and  applied  to  the  thermal 

(45  46) 
decomposition  of  Fluoroform  by  Tschuikow-Roux.    '     The  use  of  this 

procedure  is  described  for  first  order  and  an  extension  to  bimolecular 

reactions  is  straight  forward.   The  main  steps  of  this  procedure  are  now 

outlined. 

Since  the  expansion  wave  in  the  driver  section  is  isentropic  and 

reversibly  adiabatic,  then 

Y/(Y-D 
P(t)  _  fT(t) 


P5  ■  ■  ■ 


where  y  is  the  specific  heat  ratio.  Differentiating  Eq.  (D.l)  with 
respect  to  time  and  setting  T(t)  ■  T  »  the  initial  cooling  rate  is 


obtained  as  a  negative  quantity 


i]  •  f¥)  >,  H 


The  rate  of  change  of  the  pressure  (dP/dt)  is  assumed  to  be  a  constant 
thus 


P(t)  =  P5  +  (dP/dtXt-tj)  (D.3) 


where  't'  is  time.   By  setting 


4T  -  T5  -  T(t)  (D.4) 
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and  ££«  =£  (t-t.)  (D.5) 

5    5 


which  are  obtained  from  Eqs.  (D.l)  through  (D.3),  and  inserting  them 
into  the  equation  for  the  rate  of  disappearance  of  reactant  'N'  during 
the  cooling  process 


~d^t)  -  N(t)A  exp(-E/RT(t))  .  (D.6) 


the  relative  loss  of  reactant  during  the  cooling  process  can  be  written 


3 


J,n  |-i  -  -K,  E  ,  (D.7) 


where  K  =  A  exp(-E/RT  )  ,  (D.8) 

5 

N  is  the  number  of  moles  at  the  start  of  the  cooling  process  (at  tj), 
and  N  is  the  number  of  moles  at  time  't  '  when  the  temperature  has 
decreased  to  a  specific  value.   In  Eq.  (D.9)  Tschuikow-Roux  defines  AT 
at  t  ■  t,  from  Eqs.  (D.4)  and  (D.5)  as 

2 
2  T 

at--, =_ — p  (D.io) 

2T„  • 


5   RinlO 
For  a  first-order  reaction  we  have  that 
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Jta  (Hj/Nj)  =  -K(,t1  (D.ll) 

and  £n(N2/No)  =  Kun  e,  (D.12) 

where  the  subscripts  'c'  and  'un'  mean  corrected  and  uncorrected  rate 
constants.   Also,  from  logarithmic  properties 

Jin  (N,/N  )  =  Zn(N,/N  )  +  S,n  (N./NJ 
1   o         2   0  12 

and  from  Eq.  (D.7)  have 

JinO^/iy  =  in  (N2/Nq)  -  K,  e  .  (D.13) 

Combining  Eqs.  (D.ll)  through  (D.13)  the  cooling  rate  correction  is 
obtained  as 

Kc  =  (1  -  e/t^"1  Kun  (D.H) 

Determination  of  Overall  Reaction  Rate  Constant 

The  procedure  to  determine  the  reaction  rate  constant  was  performed 
assuming  first  order  behavior  found  from  reaction  time  studies.   From 
the  experimental  data  an  initial  least  squares  fit  was  performed  on  the 
variables  1/T  vs  log  K.  where 


108  K  =  log  A  "  2303RT  <D"15) 


and  K  =  Ae"E/,RT  (D.161 


E  =  energy  of  activation, 
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A  =  preexponential  term, 

T  =  temperature, 

R  =  universal  gas  constant . 

This  least  squares  produced  an  energy  of  activation  and  an 
uncorrected  rate  constant  K  .   By  applying  Eq.  (D.14)  a  new  and 
corrected  rate  constant  K  was  computed.   This  new  K  was  compared  to 
the  old  value  K   and  a  few  iterations  yielded  the  computed  value  of 

K  =■  2.5  x  10  exp(-107  (kJ)/RT)  sec-1 

2 
with  a  correlation  coefficient  of  r  =  0.8636  for  the  temperature  range 

1100-1700  K. 

These  calculations  were  computed  by  a  Basic  language  program  in  a 

HP9816  computer  which  is  shown  in  Fig.  (D.l). 
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Fig.  D.l.   Computer  Program  Listing  of  Cooling  Rate  Calculations. 


3  !         THIS  PH06BSH  CALCULATES  THE  REACTION  COOLINS  SATE 

4  !      CORRECTION  FOR  A  SET  OF  GIVEN  DATA  USINB  THE  PROCEDURE 

5  !  DEVELOPED  BV  TSCHUIKON-HOUI 


20 


OPTION  BASE  I 

DID  T3llMI,MtltllOOI ,1111001  ,KIIMI, Until (1001, !n«trllOOI,l.Dfoll30l 
21       DIN  Et(100>,P5(I00l,Dpilt(100>,Ratio(100) 

23  DIN  Initial (1001, Findl (1001, Entrqr(100>,PrttillOOI,Rco»tf  1100], Corr (100) 

30  INPUT  'INPUT  YOUR  HAIIHUH  CORR.  COEFF.  «2',Eps 

31  JasfZ./3. 

40  M.31434 

41  INPUT  'INPUT  DATAFILE  NANE  AND  TOTAL  NUB8ER  OF  DATA  PTS',Hatet,K 
50  "ASS  STORAGE  IS  ,:HP82?01, 700,0' 

60  ALLOCATE  0atal[255! 

70  Fortac:        [BABE  2I,D.5D,2I,4D,2I,D.3DE,I,4D,2I,D.3JE,2I,0.3DE,2I,D.3DE 

SO  ASSIGN  PPath.l  TO  lii>l 

90  ASSIGN  tPatti  2  TO  BUFFER  12541 

100  FOR  1=1  TO  N 

110  ENTER  *Path.l|0ato« 

120  OUTPUT  M>atii.2;Dacot 

130  ENTER  »Patn.2  USING  74>;Dtiell  1 1)  ,TSI1 1  .P3III  ,K(II,  Bpot  I II ,  Initial  II)  .Final  ill 

143  NEXT  I 

150  ASSIGN  tPath  I  TO  I 

160  ASSIGN  iPath.2  TO  I 

164  Prnl=2 

148  FOR  1-1  TO  N 

175  H(I>=-.4IT5II]IDpat(ll/P5II) 

174  ln»tflII=l/T5II! 

177  NEIT  I 

182  H«l 

183  J=0 

184  PHUT  -IOUH  CORRELATION  COEFF.  IMS' 

185  PRINT  Eps 
184  PRINT 

1B7  PRINTER  IS  1 

188  PHINI  '   R2     ITERATION  NO.' 

189  FOR  I'l  TO  N 

193  LtHolll-LGTIAISlKIIIII 

191  NEIT  I 

192  CALL  Lsquarts(Inver(ll,Lpqp(lt,N,A,8,R2) 
195  Ea=A8S(BIIRIL06(10] 

194  Entrqy[HI=Ea 
197  Pret»(HI«10-A 
191  Rcqt))UI=P,2 

199  IF  R2<Rr2  THEN  GOTO  250 

200  FOR  1=1  TO  N 

201  Dtltat(II'(2IT5(I|-21/l(Ea/IL0G(I0IIRI)tl2IT5(I)ll 

202  Et(II>(R7EalllT5(i|-2/H(I))l(l-IEIP(-EaHtItat(II/IT5II)-2IR]lll 

203  IF  Oriler=l  THEN 

204  Dtnot>I-IEt(II/0«tlltl)l 

210  KID'tlll/Ienot 

211  6010  220 


100 


212  END  IF 

21!  CALL  Secuit<Fold,Ordir,KI!>,»»elllII,initiillII,E«(I)l 

214  Kdl-Fold 

220  NEII  [ 

230  8r2=B2 

232  PRINT  82,H 

233  H=*H 

2*0  IF  S;.  EPS  THEN  6010  189 

230  FDR  I'l  ID  H 

240  li>tr<I>>[ntrl[)/.0Ml 

270  HEIT  I 

271  PRINTE8  IS  1 

273  IIIFUI  "DO  YOU  MB  10  TRANSFER  DATA  INTO  ASCII  FILES  ?  IV'l/IM)',Prb 

274  IF  Prb«0  THEN  6010  290 

200  CAa  "ransi  Inver  i'l i .  Laga I II  ,H> 

290  PRINT 

300  S=I0'A 

Ml  PRINTER  IS  I 

310  PRINT  '  EiU/Mll        PREE1P.  R2         RATE  ORDER' 

320  PRINT  En<rgy(H-I],PrM«IH-ll,8cMfHH-l]TOnler 

321  PRINTER  IS  1 
330  END 

340  I!!!!!!!!!!!!!!!!!!!!!!!!*!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!1!!!!!!!! 

341  ! 

350       !  SUBROUTINES 

340       ! 

301  !!!!!!!!!•!!!!!!!!!!!!!•!!!•!!!!!!!•!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!)!!!!!!! 
370  SOS  Lnuiresdlll, (ill, N,A, 1.821 

380  Uno"0 

390  Doi'O 

400  Tres-0 

410  Cuatro«0 

420  Cinco'O 

430  laiMt 

440  FOB  1  =  1  TO  N 

450  HnilNlUm 

460  Dos'Dos*Y[II 

470  TrtS'TrntllllmiDI 

4S0  Cujtro=Cmtn>»IHII"2l 

490  NEIT  1 

500  ItMn=Uno/N 

510  Tuu'Diis/N 

520  8l=(Trei-l(lIiiolDosl/N>l/ICultri!-(llno"2/N)l 

530  FDR  I<l  TO  N 

540  Cinca=Cinco»liIlli-<«in)'2.' 

550  Seii»Seis<(imi-y«ejnl-2l 

540  NEIT  I 

570  «2"ll"2ICinco/Stii 

433  A-VMin-IBllIieanl 

434  !<II 
440  SUSEND 
450  ! 

440  '  TRANSFERENCE  OF  DATA  TO  FILES 

470  ! 

480  Sill  Tr«ni(l(l),TII),N> 
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681 

IMSS  STORAGE  IS  ':HP82?0I, 700,0' 

682 

ALLOCATE  Ditil[255J 

483 

Fonat:      HUGE  21,40. 60,21, 40. 6D,2I, 3D 

184 

INPUT  "MAT  IS  THE  NAHE  OF  YOUR  ASCII  F1LE?',N|I|1 

135 

INPUT  '00  TOT)  mi  TO  CREATE  A  NEK  FILE?  (Y=I,H>OI', Prill 

lit 

IF  Prll=0  THEN  6OT0  705 

692 

CREATE  ASCII  UaieM 

MS 

ASSIGN  ePith.l  TO  Rust 

6»4 

ASSISN  IPjtl.2  TO  BUFFER  [2541 

*H 

FOR  l"l  TO  R 

Hi 

OUTPUT  IPit».2  USING  Foriltlllll.Yill 

m 

ENTER  0Pit«.2;Dlta( 

o98 

OUTPUT  -Pat,!   kSst.l 

o99 

REIT  1 

700 

PSSI6N  *Plth_l  TO  1 

701 

ASSISN  *Pith_2  TO  1 

702 

PRINT  'YOUR  H  VALUES  HAVE  BEEN  FILED!!!' 

703 

GOTO  731 

m 

INPUT  'HOI  NANV  [ATA  PTS.  ARE  ALREADY  F1LE0?',N 

703 

ASSIGN  ePitll.l  TO  R«et 

70* 

ASSIGN  SPatk.2  TO  BUFFER  [2541 

70? 

FOR  1-1  TO  N 

701 

ENTER  tfith.liOlUI 

700 

OUTPUT  PPath  2jData( 

'10 

ENTER  SPath  2  USIRS  Foriit;llI>,Y(I> 

711 

REIT  I 

712 

ASSIGN  PPath.l  TO  t 

713 

ASSI6R  SPitll.2  TO  1 

H< 

INPUT  'HOI  RANV  REN  I-Y  VALUES  ARE  YOU  ENTERING?', H 

715 

PRINT  'INPUT  I-Y  VALUES' 

714 

FOR  1=1  TO  N 

m 

INPUT  KI*Nl,Y(IiNI 

71? 

NEIT  I 

720 

N'N*N 

721 

PURSE  Villi 

7:; 

SOTO  6)3 

751 

SUBEND 

741 

! 

751 

t 

741 

! 

:s2« 

SUB  Pro;rlA,R«> 

7854 

Rh-I 

2894 

SUBERD 

2841 

! 

7874 

!                          SECANT  KETHOD  (BOOT  FINDERI 

2884 

! 

3844 

SUI  Secint[Foll,Order,»,DMlll,InitiilI,E«l 

3874 

EpiI'.Ol 

'875 

Ao-KHIO.I/IO 

7874 

Utrtuimiu 

7.8S4 

li»0 

3894 

Ti1=Ao/Am 

7t.)4 

Al'Ao 

7914 

Ard*Tel/Ao 

3W4 

A=Aoo 
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3929  SOSUB  4100 

3430  Fold=Ses 

3932  Mi 

3933  SOSUB  4 1 00 

3934  NeMes 

3954  IF  Foid=Ne»  THEN  SOTO  4015 

3955  Tel-Ne«4Tel/INt«-Foldl 

3964  Ttl'MSIItll 

3974  IF  TeKEpil  THEN  GOTO  4015 

3984  Ii=Ii.l 

3985  PRINT  li.Tel 
3994  Fold'Hn 

4M4  IF  Ii=50  THEN  60T0  4034 

4005  60T0  3932 

4015  EdM 

4024  SOTO  4044 

4034  PRINT  'ErERATIDK  LIMIT  RESCUED,  1=50' 

4044  SOTO  40M 

4054  PB1«T  'STOP  ITERATIONS  OUE  TO  LARGE  AO-AOO  OALUE' 

4064  SOTO  4410 

4100  One=A 

4101  Uno=cInitiiir(Grder-llHD«elllllA-Kk>l(Order-t) 
4200  Tao<OSIUnoi/IEeitlOrdtr-l)l 

4300  Res=One-T«o 

4400  RETURN 

4410  SUBEND 

5873  SUB  ProqrlA.RM.Order.tt.Dnein.inituU.Eeei 

5874  3ne=-L0S(8i 

5875  I«o=10B(IIniUjlr(0rder-ll)l[«emi(A-Kl)l(0rder-l)l/iEiel(0rder-l]; 
5883  Res=OneH.06(TM> 

5903  SUBEND 
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Appendix  E:   Computer  Programs 

The  analysis  of  the  experimental  data  was  performed  with  the  use  of 
an  HP/9816  computer  in  Basic  language.   A  brief  explanation  of  the 
computer  program  follows  next  and  a  print  out  of  the  program  can  be 
found  in  Fig.  (E.l). 

Also,  the  sensitivity  analysis  discussed  in  Section  3.3.2  was 
performed  by  a  computer  program  written  in  Fortran  for  a  mainframe 
IBM/370  computer.   The  description  of  this  program  (DGEAR)  can  be  found 
in  the  IMSL  routine  handbook,  and  a  printout  of  the  program  can  be  found 
in  Fig.  (E.2). 

Reduction  of  Experimental  Data 

The  Basic  program  (Fig.  E.l)  for  data  reduction  was  divided  into 

three  cases:   1)  Perkin-Elmer  analysis  for  C.  and  C.  compounds,  2) 

Carle,  for  C  through  C  compounds,  and  3)  Tracor,  for  products  with 

molecular  weight  higher  than  toluene.   The  program  requires  the 

information  shown  in  Fig.  E.l (a)  for  Carle  and  Perkin-Elmer  data.   The 

measured  areas  from  the  gas  chromatographs  (GC)  is  then  analyzed  by  the 
program  using  the  equation 

f  C  1   ""c  Pf 

u  =  Area  x  \- — -  -—-  (v   n 

s       s    lAreal   MW  P,  l».u 

1    '    si 
c 

where  ug  is  the  concentration  of  sample,  Area  is  the  area  of  the  sample 

from  GC,  C  is  the  concentration  of  known  standard,  Area  is  area  of 
c  c 

known  standard,  P1  and  P  are  the  initial  and  final  pressures  in  the 
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test  section,  and  MW  is  the  molecular  weight  of  C  =  standard,  or  S 
sample.   Toluene  decomposition  is  calculated  from  the  equation 


nf  _  Area(s)   f 

n    Area(mix)  P  kc'-z; 


where  Area(s)  is  the  area  of  toluene  from  sample  bottle,  and  Area(mix) 
is  the  area  of  toluene  from  Ar/C  H.  mix  bottle.   For  the  case  of  Tracor 
data,  only  three  compounds  (biphenyl,  bibenzyl,  and  Diphenyl  Methane) 
were  analyzed  quantitatively.   This  analysis  was  performed  by  first 
obtaining  a  relation  between  the  compound  area  in  the  liquid 
chromatograph  as  a  function  of  the  concentration.   Once  this 
relationship  was  found  for  the  three  compounds,  the  results  were  fed 
into  the  computer  program  which  then  calculated  the  compound 
concentration  from  a  given  area  in  a  sample. 

An  example  of  the  computer  output  is  shown  in  Fig.  E.l(b)  where  the 
concentration  for  each  compound  from  the  Perkin-Elmer  and  Carle  G.C.  is 
obtained.   The  program  was  also  capable  of  storing  all  the  data  in  disc 
files  automatically  and  simultaneously  for  a  given  shock  wave.   This  was 
of  advantageous  use  when  temperature-concentration  plots  were  to  be 
performed  because  plot  programs  were  written  in  a  way  that  permitted 
retrieving  the  whole  file  from  a  magnetic  disc. 

Modeling  of  Experimental  Data 

As  it  was  mentioned  earlier,  the  modeling  of  the  experimental  data 
was  performed  by  using  an  IMSL  routine  called  DGEAR  which  uses  a 
backward  differentiation  formula,  usually  called  Gear's  stiff  method,  to 
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solve  the  set  of  simultaneous  differential  equations.   The  program 
printout  is  shown  in  Fig.  (E.2).   The  program  was  modified  so  that  in 
the  output  a  concentration-time  plot  would  be  provided  where  the 
transient  history  of  each  species  could  be  observed.   The  input 
parameters  needed  by  the  program  were:   reaction  time,  initial  species 
concentration,  number  of  steps  in  first  iteration,  accuracy  of  result 
for  each  iteration,  and  temperature.   Also,  the  program  has  two 
subroutines  (RTCNST  and  FCN)  where  the  rate  constants  and  the 
differential  equations  are  input  by  the  user. 
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Fig.  E.l.   Computer  Program  for  Analysis  of 
Experimental  Data. 

(a)  Input  Data  for  Program 

(b)  Output  Results  from  Program 

(c)  Program  Listing 


SHOCK-WAVE  DATA  SHEET 


Run  No. 

Ta«p.   15' 

K         Pf                            P.ig 

71                    Torr          ?4 

P'lK 
7ERKTN-ELHBR 

Dwell  Time                     maec          Date 

SAMPLE  AREAS 

CALIBRATION 

Range:   X 

Cyclohexane:                               X 

CfiHft    (IflAI  KPH)! 

* 

.       _                                                                Toluene  Kill                                    r 

3.58 

6.17 

CARLE 

CM  (4  .80) 

Can  Mix  34: 

C2H4(3.85) 

CH4(4.92) 

C2H2(8.10) 

C2H2(8.39) 

C3H6(9.12) 

C3H6(9.12) 

rmiwii— «cn  no) 

C3H4 (18.46) 

C3H4U7.95) 

1,3  Butadiene(20.90) 

C4H4<25.50) 

C4H2(35.70) 

Range  X 

00 
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SHOCK  NUMBER:  831 10703 


TEMPERATURE: 


CARLE  DATA 

SAMPLE 

AREAS 

X  FINAL/INITIl 

CH4(4.80) 

4.39120 

.39401 

C2H4IS.SS) 

1.13290 

.04676 

C2H2I8. 10) 

83.32280 

3.70536 

C3H6<9. 12) 

0. 00000 

0.00000 

PR0PADIENEU2. 

00) 

. 76600 

. 02356 

C3H4(17.95> 

2. 17960 

.06703 

1-3  BUTADIENE! 

20.90) 

1.SB350 

.01213 

C4H4(23.50) 

1. 16910 

. 00732 

C4H2<33.70> 

.42810 

. 00298 

PERKIN-ELMER  DATA 


PI  (ACTUAL) 


PF( ACTUAL)    TOLUENE  DECOMP.   7. INITIAL 

TOLUENE  COMP. 


INITIAL  CONC. (mol/mlts) 
7.470E-08 


FINAL  CONC. (mol/mlts) 
1.681E-08 


■/.   FINAL/INIT.  RATIO 


3.38        S.5180 
6.17       40.3220 


.3281 
22.4988 


(b) 
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20  ! 

33  '        THIS  PfillBRBH  REDUCES  THE  INPUT  MIS  F  RON 

40  !        THE  CARLE,  PERKIN-ELNER  AND  TBflCOR  GAS  CHROHATOGRAPHS 

50  !        FOR  A  GIVEN  SHOCK  NAVE  TENPERATURE 

40  ! 

61  !  THE  NAIINUN  NUflSER  OF  DATA  STORED  [S  SO  PTS. 

62  ! 

70  !l!!l!!!!t!l!!!!!!!!!l!!!!!!!!ll!!!! !!!!!!!!!!!!!!!!!!!!!!!!!! !!!!!!!!!!!! 

80  [IK  Carle(9l,C(9l,N«(9>,Cilibl9>,N»i(9>,Pertinl9),Psrea(91,T:iel9) 

10  DIM  TriCO(22l,Co(22l,Bi|ll5),DMtl5>,BiD(l>,Pct6),Fcr.{6) 

VI  DIN  I»liOI,ryyliO),ly:l4l,Truir(221,Asciii[iO,3SI 

92  SIN  Cirlll9][201,Ojtall39U5J 

93  CON  IxllOO>,VllOO>,Cour,tllOOI,ZllOOI,FllOOi,Z:llOOI,F<<IOOI 

93  DATA  CH4(4.80l,r,2H4(5.85>,C2H2<8.iO>,C3H4l9.12>,PROPADIENEIl2.00),C3H4ll7.95l,l-3  BUIAEIENE I2< 

96  READ  CirUIII 

97  DATA  PE047,PE054,PE090,PE1I6,PE145,PE168,PE207,FE358,PE6I7 

98  DATA  CH4,C2H4,C2H2,C3H4,PR0P,C3H4,BUT,C4H4,C4H2 

99  DATA  TROOO,TR104,TSI28,TRI68,TR173,I«184,T»232,TR241,TP,093,tP,i041TR124 

100  DATA  TR148,TR20i,TR24a,TR285,TR329,TR3a9,TR445,IR493,TR542,TR667 

101  READ  Ditlflll 

102  PRINT 

103  PRINT 

104  INPUT  '   INPUT  THE  SHOCK-NAVE  NUNSER,  DNELL  TINE  AND  TENPERATURE  TS'-,Nui,Dt,T5 

105  INPUT  '   INPUT  PSITORRI.AND  TINE  lisic)  IT  TAKES  TO  P3  TO  BE  KALF(T12)',P3,T12 

107  PRINT 

108  PRINT 

109  PRINT  ■      INPUT  A  NUNSER  FOR  DESIRED  ANALYSIS:' 

120  PRINT  "  2=PERKIN-ELNER" 

130  PRINT  ■  S-TFtACOR1 

140  PRINT  "  4»CARLE,PE«KIN-ELNEFf 

160  INPUT  taker 

191  DATA  0.47, 0.34,0. 90, 1.16, 1.45,1. 68, 2. 07,3.5a, 6. 17 

192  READ  TiKltl 

193  DATA  92.13,78.114, 10BI 

194  READ  N«t, Nab, Spin 

196  DATA  .2076,0.1629,0,1.2409 

197  READ  Jy:  .0 

200  DATA  15.1,13.98,13.98,11.84,13.00,13.00,16.1,16.1,16.1 

210  READ  Clll 

220  DATA  16.0430,26.0382,26.0382,42.0810,40.0653,40.0653,54.0924,54.0924.54.0924 

230  READ  Nam 

240  DATA  16. 0430,28.0542, 26. 0382,42.0810,40.0653, 40. 0653, 54. 0924, 52. 0765, 50. 0605 

250  READ  M  111 

257  DATA  0.00,1.04,1.28,1.48,1.73,1.86,2.32,2.41,0.93,1.06,1.24,1.48,2.01,2.43,2.85 

258  DATA  3.29,3.89,4.45,4.93,3.62,6.67,7.00 

260  READ  Tricirll) 

261  DATA  0.0001,0.00001,0.000001,0.0000001,0.00000001,0.000000001 

242  READ  Pclll 

243  DATA  94.3584,11.3724,1.3893,0.2076,0.0354 

264  READ  Biplll 

24S  DATA  842. 2195,84. 0192,8. 1372, 1.2514,. 1429 

244  READ  Out  III 

267  DATA  879.4521,86.3488,9.3852,1.2409 

265  READ  Siblll 

Cc) 

no 


270 

Cccc=0. 

27! 

IF  Nu«ber=2  OR  Nuiber=4  THEN  GOTO  570 

273 

IF  N«»ber=3  THEN  SOTO  1022 

274 

nil ill  II II II  nil Ilttlllllllliilitiiiiiiiin ti 

2/5 

! 

277 

!                                       CARLE  PROSFtAN 

278 

• 

279 

FOR  1=1  TO  13 

230 

PRINT 

211 

NEIT  1 

290 

PRINT  '                 STABT  B»  IMITIN6  CABLE  DATA1 

300 

POINT 

301 

PRINT 

310 

PRINT  ■      INPUT  SANPLE  AREAS  (CARLE  PRINTOUT  OROEfi,  ONE  S»  ONE)  ■ 

311 

PRINT 

312 

PRINT 

320 

FOR  1=1  TO  9 

321 

PRINT  CarlKIl 

330 

INPUT  Carlelll 

340 

PRINT  I.Carleii) 

330 

NEXT  I 

3S1 

PRINT 

332 

INPUT  '  INPUT  THE  RAN6E  FOR  SANPLE1, Rl 

3S3 

PRINT 

354 

FOR  I>1  TO  9 

355 

CarIelII=Cjrl«(I)IRl 

356 

KEIT  I 

357 

FOR  M  TO  17 

33! 

PRINT 

M 

NEIT  1 

362 

PRINT  -INPUT  CALIBRATION  AREAS  ONE  U  ONE  (CH4,C2H2,C3H6,C3H4r 

363 

PRINT 

361 

PRINT 

370 

5o«'0. 

380 

FOR  1=1  TO  4 

361 

IF  1=2  OR  1=3  THEN  EOTO  420 

353 

PRINT  Cdrlllll 

390 

INPUT  Clliblll 

4u0 

PRINT  I, Ciliblll 

410 

S«ib=Suib»(CiliblII/Clin 

420 

NEIT  I 

430 

Cllibl2l=Cilibl3> 

440 

Cllibt5l=Cjlib(i> 

450 

Cilibl7>=ICI7>ll00ISlnbl/l69.l769) 

460 

C'libl8l=Cilibl7l 

470 

Cil>bl9l=Calibl8l 

471 

PRINT 

472 

INPUT  'INPUT  RANGE  FOR  CALIBRATION-, R2 

473 

PRINT 

473 

PRINT 

474 

FOR  1=1  TO  9 

477 

CiIiblII=CiliblI)IR2 

478 

NEIT  I 

479 

FOR  1=1  TO  9 

480 

Clll'IICarlHIIICil]|H»IIII7(N«ilI)lCilib(IIII/(l00lTol«eiul 
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<81      NEXT  I 

483  FOR  1  =  1   TO  3 

(St      P8INT 

'35  KIT  1 

'34      IF  »uiber=4  THEN 

«7  PRINT  USIN6  15I,15A,8B,5I,14A,4B1;1    SHOCK  NUTtlER:  ■,Nui,-TENPE«ATURE:  -,15 

488  PRINT 

489  END  IF 

490  PRINT  ■  CARLE  DATA1 

491  PRINT 

W  rami  "  sanple  areas       I  fihal/ihitial  ratio/dneli.1 

493  PRINT 

500  FOR  W  10  9 

520     PRINT  USIN6  1IOI,20«,2I,40.5017I,2B.58-!C4rlllII,Cjrl«(II,ClII 

530  NEXT  I 

531  IF  Ans=l  AND  NcuberM  THEN  SOTO  903 

532  INPUT  -DO  TOO  NANT  A  CONPLETE  PRINTOUT?  (T=l,N>Or,A»s 
537  IF  ftos=l  AND  lkut<r>4  THEN 

338  PRINTER  IS  701 

559  8010  483 
540  END  IF 

543  IF  Nuner-4  THEN 

544  PRINTER  IS  I 

546  SOTO  1343 

547  END  IF 

548  "-!<"  """inn, mimm, ,,,,,,,,,,,,,,,,,,,,,,,,1,1,, ,,,,,,, .,„ 

550  ! 

551  !  PERKIN-EU4ER  PROGRAN 

560  ! 

570  FOR  I»l  TO  13 

580  PRINT 

591  REIT  I 

590  PRINT 

631  PRINT 

640  PRINT  -INPUT  THE  PERKIN-EU1ER  SANPLE  AREAS  FOR  8ENIENE  AND  TOLUENE1 

660  FOR  [.8  TO  9 

670  PRIKT  TiMIII 

680  INPUT  Piriilll 

681  PRINT  I.PirulI) 

690  NEIT  I 

691  PRINT 

692  INPUT  'INPUT  THE  RANSE  OF  THE  SANPLE1, R3 

693  PRINT 

694  FOR  1-1  TO  2 

695  H'I»7 

697  Plr!i(Hl=PjrH(H)IR3 

698  NEIT  1 

699  FOR  I«l  TO  20 

700  PRINT 

701  NEIT  I 

702  INPUT  -INPUT  TOLUENE  HIITURE  AND  BENZENE  CALIB.  AREAS1, Tol, Bin: 

703  INPUT  ■  INPUT  THE  RANSE  FOR  THESE  AREAS1, R4 

704  Tol«TollR4 

705  !eni-B(»!lR4 

770  Perl!in(8l=flppalP*real8)/(Beni4I00001 

773  Tolu6n«-.ll55ITol/Be»i 
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775  Perkin(ai  =  (Perkin  (81  /(Toluene!)  1100 

780  PUNT 

790  PRINT 

800  PRINT  'INPUT  THE  FINAL  SHOCK  TUBE  PRESSURE  IPSI6P 

810  PRINT  -UNI  THE  INITIAL  RESION  1  PRESSURE  (T0RR1  - 

820  INPUT  PI, PI 

830  PI-IPIII. 036251-1. 93B67 

840  Pl'IPll. 891981*45. 12204 

850  A=Pltl4.7 

860  D"PII. 01934 

863  Ed»A/B 

870  Rat=Parea(9IIEd/Tol 

B71  Dndt=IP5-PmiOOO/T12 

872  Tolinil"ToluenelPII5.3807E-IO 

873  Tollin»TolinitlRit 
87!  Kl  =  -10004L0GIRjtl/0t 
876  Et  =  Btl.001 

878  INPUT  ■  DO  YOU  BANT  TO  TRANSFER  1ATA  FOR  RIH  OROER  CALCULATIONS?  [M,M',MI 

87?  IF  Prb2-1  THEN 

880         INPUT  'PRINT  'RIORDER'  AT  THIS  TIHE'.Natel 

S9I  CALL  RorderlDt,T5,P5,Kl,lpdt,rolinit,Tollin,Na««l 

882  ENB  IF 

883  PRINT  -K1=";KI j  ■  DPDT-'tOpdt 

884  PerHn(9l"Ratll00 
883  FDR  1=1  TO  10 
a»o    PRINT 

WO  NEXT  I 

901  IF  AOS'I  AND  Nuil>er>4  THEN  PRINTER  IS  701 

903  FOR  1=1  TO  3 

904  PRINT 

905  NEIT  I 

'06  PRINT  ■  PERHN-ELNER-DATA' 

907  PRINT 

908  PRINT 

"0  PRINT"         P1IACTUALI    PFlACTUALI   TOLUENE  »ECO«F.  I1NITIAL" 

'"  p™  "  TOLUENE  COUP." 

913  PRINT 

920  PRINT  USING  ■!3I,4D.4D,4I,4D.4D,8I,40.4D,4I,SD.4D-;Pl,Pf , Rat, toluene 

921  PRINT 

'22  PRINT-  INITIAL  COHC.Ieol/eltsi  FINAL  CONC.  laol/tlts!" 

923  PRINT  USIN8  "tai,D.3DE,19I,D.3uE";Tolinit., Tallin 

930  PRINT 

'50  PRINT  ■  TIltE     AREA    1  FINAL/INK.  RATIO/DNELL' 

951  PRINT 

960  FOR  [«i  TO  9 

970     PRINT  USING  '131,40.21,41, 4D.4D,4I,4D.40"lTiK(II,Parea(II,PerkinlIl 

980  NEIT  I 

981  IF  Cccc=3  THEN 

982  PRINTER  IS  1 

983  PRINT  'END  OF  PROGRAH  REACHED" 
'85         GOTO  1419 

'86         END  IF 

"3  IF  Nuiber'2  THEN 

''4  INPUT  'DO  YOU  KANT  TO  HAVE  »  PRINTOUT?  IM/MIMn 

"5  IF  An'.l  THEN 

"6  PRINTER  IS  701 
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997 

Cccc=3 

9)8 

SOTO  903 

999 

END  IF 

ION 

eno  IF 

IMI 

[F  Ansa)  AND  Nu«4er"2  THEN 

IMS 

PRIMT  -EM  OF  PR06RAH' 

1403 

SOTO  1419 

1004 

EW>  IF 

lOOo 

IF  Ans.l  AND  NuiterM  THEN 

1007 

PRINTER  IS  1 

1009 

SOTO  1343 

1010 

END  IF 

1011 

IF  NoiberM  THEN 

1012 

PRINT  'PUSH  CONTINUE  TO  RESUHE' 

1014 

PAUSE 

1015 

SOTO  179 

1014 

END  IF 

1013 

( nt i iiiiiiiiiiiiitiiiin 

1019 

! 

:o;o 

!         TRACOR  PACSRAN 

1021  ! 

1022  PRINT 

1044  FOR  1=1  TO  19 

1050  PRINT 

1051  REIT  I 

1052  IF  tiu.oer=5  THEN  SOTO  1056 

1054  INPUT  -INPUT  PRESSURE  PIIT0RRI-.P1 

1055  Pl*<PH.!9l9flH45.12204 

1056  PRINT 

1060  PRINT  "INPUT  SANPLE  AREAS  ONE  BY  ONE1 

1061  PRINT 

1070  PRINI  'START  KITH  T-150  C" 

1071  PRINT 

1000  FOR  I-l  TO  21 

1001  IF  M  THEN 

1002  TranIII=0. 

1053  SOTO  1100 
1004         ENO  IF 
1096    PRINT  Tncarlll 

1090  INPUT  Trjcolll 

1091  IF  H  THEN  PRINT  'INPUT  AREAS  FOR  W20  C" 

1100  NEIT  I 

1101  PRINT 

1102  PRINT 

1103  INPUT  'INPUT  THE  RANSE  OF  THESE  AREAS1, Rr5 
1108  PRINT 

1110  PRiNT  'INPUT  CALI!  AREAS  FOR  BIPHENTL,D-HETHANE,BI&ENZtL- 

1111  PRINT  -NOTE:  THE  SAN6E  SHOULD  BE  THE  SANE  AS  FOR  SAMPLE  AREAS1 
1120  INPUT  Biph,I»tt,liDn 

1112  PRINT 

1152  INPUT  -INPUT  THE  TOLUENE  CALIBRATION  AREA  AT  150  C  ANt  ITS  SANGE'.Tratol.firs 

1162  BirU'PIITritollRrtl. 06893 

1163  PRINT  'TOLUENE  INITIAL  AREA' 
1112  PRINI  Sim 

1190  Ai'Iy:lll-!:pn 

1191  FOR  I«l  TO  5 
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1192  V«tll)=lip(ll 

1193  NEIT  1 
UN    MyjUI 

1200  FOR  1=16  10  It 

1201  IF  1  =  13  THEN  SOTO  1217 

1204  CULL  LoqilAi,I,VKtlll,Pclll,Pcclll,lll6,II 

1205  PRINT  •         A  !■ 

1206  PRINT  lt,S 
1208    IF  11=0  THEN 

120?         PRINT  ■  DEVIATION  IS  LARGER  THAN  !02  FOR  ADJUSTING  TRACOR  CURVES" 

1210  PRINT 

1211  PRINT  USINS  ,5I«,2J.2JI,4»'|"  PR06RAN  HAS  CONFUTING  CONCENTRATION  FOR  TiHE  i 

1212  GOTO  1217 

1213  END  IF 

12H  IF  Trico(I)=0  OR  Tritolll»I  THEN  GOTO  1217 

1215  PoflLOGITricolIH-NI/S 

1211  CoillMO'lPntl 

1217  H=[-14 

1222  IF  1=19  THEN  GOTO  123G 

1223  MyzlHI 

122!  IF  1=16  THEN  A.=Iyz(H)-OutH 

1224  IF  1=18  THEN  Ai=IyzlH>-iil>en 

1228  IF  1  =  16  THEN 

1229  FOR  !•!  TO  5 

1230  VecttJI'IietlJ! 

1231  SET!  J 

1232  END  IF 
12)3    IF  1=18  THEN 

1234  FOR  J=l  TO  4 

1233  VectlJ)>liblJI 

1236  NEIT  J 

1237  END  IF 

1238  NEIT  I 

123)    FOR  1=1  TO  21 

1240  IF  l"3  THEN  GOTO  1243 

1242  Tr.MlI>=TricoimRr5/Bjre. 

1243  NEIT  I 

1244  F8R  1=1  TO  20 

1245  PRINT 

1234  NEIT  I 

1255    IF  Ans=l  AND  Nmber=5  THEN  PRINTER  IS  701 

1257  FOR  M  TO  3 

1258  PRINT 

1259  NEIT  I 

1260  PRINT  ■  TRACOR  DATA' 

1261  PRINT 

1264    PRINT-  TINE       AREA       CONC.  iml/il)  --  TINE         AREA  CMC.  Iiol/ill- 

1274  PRINT 

1275  PRINT  •  M50  p 

1277  Nn=l 

1278  »i=2 

1279  FOR  Ml  TO  Nl 

1282  PRINT  USING  ■6I,DD.DD,3I,D.3DE,3I,D.3DE,3I,r!Tr.cjrIII1trscoiIl,CiiiI] 

1283  NEIT  I 
I2B4  PRINT 
12G5    Nn=fliM 
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1281  IF  Nn=9  IHEN  PRINT  '  1*220  C 

1288  *I1tl*2 

1289  IF  «i-23  OR  Nl<23  THEN  GOTO  1279 

1290  IF  Nwtir-3  THEN  >it-Rit«l  H4.5071ITriCOI3H8r5l,  ITritollRrUF  1. 1 

1291  IF  Cccc=4  THEN 

1292  PRINTER  IS  I 

1293  PRINT  'END  OF  PROGRAN' 
1295        GOTO  1417 

1294  END  IF 

1297  PRINT 

1298  IF  A»s=l  AND  Nlllbir'S  THEN  GOTO  1318 

1299  IF  Nuaber'5  THEN  GOTO  1321 

1301  IF  NutttrM  THEN 

1302  PRINTER  15  1 

1303  SOTO  1419 

1304  END  IF 

13D5     INPOT  "DO  YOU  NANT  TO  HAVE  A  PRINTOUT?  IV-l/N-O)  ■  ,(U1I 
1301    IF  Ani=l  AND  Nuibir-3  THEN 

1307  PRINTER  IS  701 

1308  CcoM 

1309  SOTO  1257 

1310  END  IF 

1311  IF  Ans=l  AND  (Nuitefl  OR  NuiDer=4J  THEN 

1312  PRINTER  IS  701 

1313  SOTO  490 

1314  END  IF 
1313    PRINT 

1311    PRINT 

13ia  IF  Ans-I  AND  Nuiker=5  THEN  PRINT  OSINS  •30MD.5DV       TOLUENE  DECOHP.  CORRECTEE  • 

1319  PRINTER  IS  I 

1320  IF  Nmbtr-5  THEN  GOTO  1343 

1322  PRINT  'END  OF  PR06RAN  REACHED' 

1323  ! 

1324  PRINT 

1325  IF  Ans'l  THEN  SOTO  1410 

1321  INPUT  "DO  YOU  NANT  TO  HAVE  A  COHPLETE  PRINTOUT?  H'l/N=0l',Ao5 

1327  IF  Ans'l  THEN 

1328  PRINT  '     LINE  OP  RIBBON  TO  PAPER,  AND  PRESS  CONTINUE  KhEN  READV 

1329  PAUSE 

1330  PRINTER  IS  701 

1331  IF  N«i«f  =  l  THEN  SOTO  483 

1332  IF  NuDer*2  THEN  SOTO  903 

1333  IF  NuiDer-3  THEN  GOTO  1257 

1334  IF  NiiMirM  THEN  GOTO  483 

1335  IF  NuiDn-5  THEN  SOTO  483 
1331    END  IF 

1337  PRINI 

1338  M!l!!!!!!!!!!!!!!!!!!!i!!!>!!!!!!!!!!!!!!!!l!l!!!!!!!!!!!!!!!!!J!!<!!!!!! 

1340  ! 

1341  !  PREPARIN6  DATA  TRANSFERENCE  TO  ASCII  FILES 

1342  '. 

1343  INPUT  'DO  YOU  NANT  TO  TRANSFER  DATA  INTO  ASCII  F IL£S?(V= 1 /N-OI -. Prb  1 

1344  IF  PrDM  THEN  BOTD  1419 

1345  IF  Prllil  THEN  INPUT  'ARE  YOU  CREATIN6  HEN  ASCII  FILES?<Y'l/N=OI',Prt 
1341  IF  Prll-I  AND  Prb=0  THEN 

1347  INPOT  'HON  HANI  DATA  PT5.  ARE  ALREADY  IN  FILE?', Chris 
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1348 

H=Chris*i 

1349 

END  IF 

1350 

IF  PrDl'l  AND  PrtFl  THEN  K'l 

1351 

! 

1352 

!                          PREPARING  PERKIN-ELHE8  DATA 

1353 

! 

1354 

FIR  1-1  10  9 

1355 

IF  J=8  08  J=9  THEN 

1354 

AsciiilN,J>=Perlin.U) 

1357 

GOTO  1360 

1358 

END  IF 

1359 

AiciiilN,J)»Pirei.(J> 

1360 

NE1T  J 

1361 

1362 

!                           PREPARING  CARLE  DATA 

1343 

1314 

FOR  ]>I0  TO  18 

1345 

Hh-J-9 

1366 

AlciiilN,J>=CiHM 

1347 

NE1T  J 

1348 

1369 

!                       PREPA8ING  TRACOR  DATA 

1370 

1371 

FOR  J»H  TO  39 

1372 

IF  J-19  THEN 

1373 

AsciiilN.J)*). 

1374 

GOTO  1383 

1375 

END  IF 

1374 

IF  1=34  OR  J>35  08  J*37  THEN 

1377 

HIW-18 

1378 

AsnulN,J]=CotHM 

1379 

GOTO  1383 

1399 

END  IF 

1381 

Jj>J-IB 

1382 

Asciii(N,JI»Traco(Jjf 

1383 

NEIT  J 

1384 

! 

13B5 

!                      REARRANGING  DATA  TO  BE  FILED 

1384 

! 

1387 

IF  Prbl-l  AND  Frb-1  THEN 

1386 

FOR  J=8  TO  IB 

1391 

Nue('DatillJ) 

1592 

In(ll>T5 

1393 

Tyylll'AsciiilN.JI 

1394 

'orro=l 

1394 

CALL  TrinsUu(t),Yyy(t),Naief,PrM,Zorro} 

1397 

NEIT  j 

1398 

END  IF 

1399 

IF  Prbl«l  AND  Pr»=0  THEN 

1400 

FOR  J=B  TO  18 

1403 

imflitam 

1404 

IKIN1-T5 

1405 

Yyvim=«scinlN,J) 

1407 

CALL  Trinslln  III, Iyy<I>.N»e',Prb, Chris, N> 

1408 

NEIT  J 

1409 

'ND  IF 
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1410  PRINT 

1411  IF  Prbl=0  THEII  8010  1419 

1412  !!!!!l!!!!!!!!!!!i!!!!!!!!!!!!!l!!!!M! !!!!!!!!!! III! !!!!!!! 

1413  ! 

1414  !  SUBROUTINES 

1415  ! 

1416  I!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!1!! ' Him 

1417  ! 

1418  !       AUIILIAR  ROUTINE 

1419  PRINT  "END  OF  PR06RAN1 

1420  END 

1421  SUB  TrKOf(«l,VKt(ll,Pcll),Pcc(ll,ll,G,Kid) 

1422  FOR  1=1  TO  Nad 

1423  IF  «i>()  THEN 

1424  V«r=VsctlII-Ai 

1425  IF  Vir<0  THEN 

1421  PRINT  'DIFFERENCE  BET1EEN  CALIB  AREA  IN 

1427  PRINT  -IN  THE  PR06RAN  IS  TOO  LARGE  ' 

1428  PRINT  ■CONPUTE  CONCENTRATION  It  HAND1 

1429  N=0 

1430  8=0 

1431  60TD  1444 

1432  ENB  IF 

1433  »ectlI)=LOSIVir> 

1434  END  IF 

1435  IF  AKO  THEN  VectlIl=L06IVect(IltABS(Aj]l 
1434          PcclIllSKPtclIll 

1437    NEIT  I 

1441    CALL  LsgiinslPcclll,Vectlll,ltiil,N,SI 

1444    SUBEND 

1443  ! 

1444  !  PERFORltS  LEAST  StUARES 
1*47    I 

1448  SUB  L'mureslllll.mi.H.A.B! 

1449  Uno=0 

1450  [<B=0 

1451  lrti-1 

1452  :uitr.i=,| 

1453  Cinto-O 
1461  3«is=0 

1471  FOR  W  TO  * 

1481         Uno=UiM»III) 

1491  hi'DlKKI) 

1501         Tres'TreiKIlIIITIDI 

1511         Cuatro=CiiitrGtUII)A2) 

1521  NEIT  I 

1531  Iiun'Uno/N 

1541  Tiein-Oos/N 

1551  Bl=>  (Tres-4  (UnolOosI  /N>  1  /  ICuitro-  (Uno'2/«l ) 

1541  FOR  1=1  TO  N 

1371         Cinco=Cincotl(<IIl-lMjn)-2) 

15S1         Seis=Seis*t(III)-YtMn)"2] 

1591  NEIT  I 

160!  R2=Br2ICinco/Stii 

1611  IF  ABSIR2K.7  TNEN  OOTO  1644 

1621  FOR  1=1  TO  N 
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Ii3l        rylll-ViuaKliKKll-Iiianl) 

1441  NEIT  I 

1442  GOTO  1664 
1444  PRINT 

1451  PRINT  -VALUE  112  IS  TOO  U»;  THIS  LEAST  SQUARES  TECHNIOUE  IS  HI  APPROBATE" 

1441  PRINT  ■   82  ■ 

1662  PSIHT  R2 
146]  PRINT 

1664  A'TMan-ltlllalanl 

1665  8=81 
1670  SHEENS 
1471  ! 

1672  !         TRANSFERS  DATA  INTO  FILES 

1673  ! 

1675  SOB  TranslInill.yyylU.NaMt.PrM.NI 

1676  NASS  STORAGE  IS  ':HPG2?0I, 700,0" 

1677  ALLOCATE  Ddtoli2561 

1679  Foriat:  IKAGE  21,60.60, 21, 6D. 60 

1681  IF  Prb=0  THEN  SOIO  1702 

16E7  CREATE  ASCII  Naut.N 

1638  ASSIGN  JPitll.I  TO  Hurt 

166?  ASSIGN  ePith.2  TO  BUFFER  [2561 

1690  FOR  1=1  TO  N 

1692  OUTPUT  fPjth_2  USIN6  FafMt!lxx(I>,fyy(II 

1693  ENTER  SPath  2;0ato« 
1494  OUTPUT  Watii  l;Datoi 
1695  NEIT  I 

1677  ASSIGN  ?Pdth  1  TO  I 

169G  ASSIGN  »Pdth  2  TO  I 

Um  PRINT  USING  ■S(M,I,IM-l"rOUR  H  WLUES  HAVE  BEEN  TRANSFERRED  TO  FILE;",Naiajl 

1700  GOTO  1722 

1702  ASSIGN  ePdth.l  TO  Nliet 

1703  ASSIGN  SPatti  2  TO  BUFFER  [2561 

1704  FOR  1=1  TO  N 

1705  ENTER  BPath  l;Datol 

1706  OUTPUT  !Patt.2;0atot 

1707  ENTER  Wath  2  USING  ForMt;Iii(II,Vyyl[l 
17CB  NEIT  I 

1709  ASSIGN  Watll  1  TO  I 

1710  ASSIGN  »Pdth  2  TO  I 

1711  PURGE  Naiel 
1720  GOTO  1697 
1722  SUBENO 
17J2  ! 

1742  !  AIIIILIAR  ROUTINE 

1752  ! 

1762  SUB  Loop(Aa,I,llKtin,PcllT,Pcclll.lt,G,Il 

1772  Nad=5 

1773  IF  Ml  THEN  NadM 
1775  FOR  Hh»l  TO  5 

1774  IF  1*14  THEN  Pec IHhl=Pc IHh*l ) 

1777  IF  |.|7  THEN  PcclHhl-PdHlil 

1778  IF  I.|9  AND  Hh=5  THEN  GOTO  1780 

1779  IF  1=19  THEN  PcclHNI'PclWll 

1780  NEIT  Hlt 

1792  T«stl=ABS(Aa/IllT0O 
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1742    SUB  Loii|ilAa,I,Vectll),Pcll),Pcclll,N,G,l> 

1794     IF  7«sti>5.0  THEN  CALL  Tricorlflj, Vect ( t) ,Pc (l(,Pcc (II ,M,6,»»li] 

1797  IF  Teitl<5.0  AND  1=14  THEN 

1798  "=14.4380 
179?                           6=1.9840 

1800  END  IF 

1801  IF  TestKS.O  AND  [=1?  THE* 

1802  »=1S.47354 

1803  8=2.19098 

1804  END  IF 

1805  IF  TestKS.O  AND  1=17  THEN 
1804  »=15.ll?3 

1807  8=2.1308 

1808  END  IF 
1810    SUBEND 

1820    ' 

1830    ! 

1840    ! 

1B50    SUB  Rord-rlUno,Doi,Trei,Ciiatrci,Cinco, Sets, State, Naiel) 

1840    CON  I<(l), Till, Countlll.illl, Fill, 7:111, Ff(l) 

1870  NASS  STORAGE  IS  "1HP82901, 700,0' 

1871  ALLOCATE  Datai[255] 

1B72  Foriat:    INASE  2I,D.5t,2I,40,2I,D.30E,I, 40,21, 1.30E,2I,D.3DE,2I,D.3DE 

1873    Cmintei-1 

1890    INPUT  -DO  »UU  NANT  TO  CREATE  A  HEN  FILE?  (V=l,N=or,Prlt 

1900    IF  Pra=0  THEN  SOTO  2040 

1911    N»l 

1920  PRINT  -INPUT  lWLL(iMtl,T3', P5(I«Tl,l(l,NP/dT,[I]c,[I]f 

1930  FOR  1»1  TO  N 

1941  I.(I)=Uno 

1942  r(I)=0M 

1943  Cotintlll-tres 

1944  KIKutra 
1944          FULCiitco 

1947  Z:l!l=Seis 

1948  Ffll!=Siete 

1950  PRINT  USIN6  Foriitjlulll, Til), butt II I, ill), Fill, IllII.BIII 

1940  NEIT  I 

1970  CREATE  ASCII  Naael,N 

1980  ASS1SN  "Patlt.1  TO  Null 

1981  ASSIGN  BPatlt  2  TO  BUFFEF.  [254] 
1990  FOR  1=1  TO  N 

IHIKI  OUTPUT  IPatl.2  USING  ForMt  j  I<  (I> ,  If  (1 1  .Count  ( I) ,  I  (I  I ,  Fill,  1:11),  F(II) 

2001  ENTER  SPatlt^iData* 

2002  OUTPUT  SPatli  l;Datal 
2010  NEIT  I 

2020  ASSIGN  SPatll  1  TO  I 

2021  ASSIGN  BPath  2  TO  I 

2023  PRINT  •YOUR  M  VALUES  HAVE  BEEN  FILED"" 
2030  GOTO  2190 

2040  INPUT  'HON  HANK  DATA  PIS.  ARE  ALSEADT  FILED'', N 

2041  N-l 

2050  ASSIGN  >Patli.l  TO  Haiet 

2051  ASS16N  SPatlt  2  TO  BUFFER  [254] 
204O    FOR  |«|  TO  N 
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SMI  ENTEN  "Pith  KSitil 

20«2  OUTPUT  SPilft  2;0itat 

!»'■«  am  ?Plth_2  USiN6  ForMt|IitII,m>,CoiintlI>, III!, Fill, ZKD.Fflll 

::bo  seit  i 

:i)6l  ASSieil  PPltn.l  to  I 

2062  RS5I6N  With  2  '0  I 

20M  FOB  1=1  10  »" 

20f5  J=N»I 

2*fl  Ix(J)=Uno 

2098    (Ui=jM 

2C?9  Count  UWrn 

2100  IIJKuitro 

2!0I  FUKinco 

2102  ZzlJI-Stii 

2104  KMi*Hit» 

210?  PRINT  USING  Forut|Ii(JI ,TIJ).Cw»t(JI, III), FIJI, IllJI, FUJI 

2108  NEII  1 

2109  PRINT 

2110  PRINT 

2111  «=IHN 

2112  PRINTER  IS  1 

2113  FOR  1. 1  TO  N 

2114  PHUT  USING  Foraitlb U),T([] ,Cannt(l>,ZIIl,F(I) ,Ii(I).F(III 

2115  NEIT  I  ' 
2114  PRINTED  IS  1 

2117  PURSE  Duel 

21  IS  SOTO  H70 

2K0  3USEN0 
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Fig.  E.2.   Computer  Listing  of  DGEAR  Modified  Program. 


file:  toluene  Fortran  a   Kansas  state  university  vm/sp  cms 

INTEGER  IY. METHf MITER , INDEX , I WK( 50 ),IER,K, INT 
REAL  T,EPS,TMAX,H,TCALL 
EXTERNAL  FCN.FCNJ 
DIMENSION  Y(20),WK(500> 
COMMON  DK(40> 

CHARACTERS  ITITLE (76)  r  CHAR (20)  /'A'  ,  'B't'C'l  'D'  t'E't'F'  ,  *B'»  'H'  I 
1   'I'.'J'.'K'r-L'.'M'.'N'.'O't'p'F'a'.'R'r'S-.'T'/ 
SO     READ  (5.1)  IY, IK, ITITLE 
IF  (IY.EO.O)  STOP 
READ  (5f3>  TEMP 
READ  (5f2)  MITERrlNDEX 
READ  (5.-1)  (Y(I),I=1,IY> 
READ  (Sfil)  METH, INT.TMAX.EPS 
T=0.0 

H=TMAX/INT 
CALL  RTCNST(TEMP) 
WRITE  (6,5)  ITITLE 
WRITE  (6,22) 

WRITE  (6,21)  ( I,DK(I ) ,1=1, IK) 
IF  (METH.EO.l)  THEN 

WRITE(6,30> 

WRITE(4,27> 

WRITE(4,29>  H 
ELSE 

WRITE(6,30> 

WRITE(6,28> 

WRITE(6,29>  H 
END  IF  . 

WRITE(4,32)  MITER 
WRITE  (4,31)  TMAX,EPS 
CMAX=Y(1> 
DO  500  1=2, IY 
300    CMAX=AMAX1(CMAX,Y<I>> 
AN0RM=1.0/CMAX 
CALL  GRAPH1  (0.0,  TMAX,  0.0, 1 .0,  ITITLE, '  ') 
DO  699  1=1, IY 

699  CALL  GRAPH2  (0.0,ANORM*Y(I ) ,CHAR( I > ) 
WRITE  (6,24)  T 

WRITE  (4,21)  (J,Y(J),J=1,IY) 

INT=0 

TCAL=T 

700  INT=INT+1 
TCAL=TCAL+H 

CALL  DGEAR(IY,FCN,FCNJ,T,H,Y,TCAL, EPS, METH, MITER.  INDEX,  IWN,WK,IER> 
IF(IER.GT.12S>  GO  TO  9999 
DO  9898  11=1, IY 
9898      Y(II)=(Y(II)*100/2.87E-07) 
WRITE  (9,24)  T 
WRITE  (9,21)  (J,Y(J),J=1,IY) 
WRITE  (6,24)  T 
WRITE  (4,21)  (J,Y(J),J=1,IY> 
IF(INT.GE.IOO.OR.T.GT.TMAX)  GO  TO  2000 
DO  1000  J«1»IY 
1000   CALL  GRAPH2  (T,ANORM*Y( J) ,CHAR( J) > 
GO  TO  700 
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file:  toluene  Fortran  a   Kansas  state  university  vm/sp  cms 

2000  write  (6.34)  int 
call  graph3 

GO  TO  50 
9999   CONTINUE 

WRITE  (6,33) 

WRITE  (6, 25>IER,METH, MITER, INDEX 

WRITE  (6,23>TCALL,H 

WRITE  (6,21)  (J,Y(J),J=1,IY> 

STOP 

1  FORMAT  (2I2,76A1) 

2  FORMAT  (212) 

3  FORMAT  (1E10.3) 

4  FORMAT  (3E10.3) 

5  FORMAT  (1H1,30X, 'KINETIC  MODELING  PROGRAM' ,///,20X,76Al ,///> 
11     FORMAT  (2IS,2E10.2> 

21  FORMAT  (5(I4,1PE1,4.S>) 

22  FORMAT  (10X,'RATE  CONSTANTS',/) 

23  F0RMAT(/,10X,'TCALL=',1PE12.5,'SEC',4X,'LAST  H=' .1PE12.5, 'SFC ) 

24  FORMAT  (//,21X,'AT  TIME=' ,lPE13.5,/> 

25  F0RMAT(//,10X,'IER=',I3,'   METH=',I3,'   MITER*', 13,'   INDEX*', 13) 

27  FORMAT  (//,10X, 'VARIABLE-ORDER  ADAMS  PREDICTOR  CORRECTOR  METHOD') 

28  FORMAT  (//, 10X, 'STIFF  METHODS  OF  GEAR  -  DIFFERENCE  FORMULAE') 

29  FORMAT  (/,10X, 'FIRST  STEP  SIZE  '.1PE12.5,-  SEC) 

30  FORMAT  (/,10X, 'VARIABLE  STEP  SIZE  CALCULATION') 

31  FORMAT  (//,10X,'TMAX=',1PE12.5,'     EPS=' , 1PE12.5, '   RELATIVE  ERRO 
1R  OF  BASIC  METHOD',///) 

32  FORMAT  (//,10X, ' ITERATION  SCHEME  *',I2,'  APPLIED') 

33  FORMAT  (//, 10X ,' TERMINAL  ERROR  OCCURED-  HELPFUL  VALUES  LISTED  BELO 
1W',//,10X, 'PLEASE  CHECK  ERROR  PARAMETER  CIER3  FOR  PROBABLE  CAUSE') 

34  FORMAT  (//,10X, 'NUMBER  OF  DATA  POINTS  IN  TIME  GRID=' ,13,///) 

SUBROUTINE  RTCNST(T) 

COMMON  DK(60) 

Rl=8. 31434*0. 001 

R2=l. 98545*0. 001 

DK(l)=2.7*10**14.0*EXP(-73.4/(R2*T>) 
C     BK(2>=10**(9.3-(1.5/(2.3*R2*T*1000>>) 

DK(2)=0. 4*10**13. 0*EXP(-3.7/(R2*T>) 

DK(3)=0.0000760*EXP(-340.0/(R2*1000*T))*(T**5.5) 

DK( 4>=8. 3*10**17. 0*EXP( -300. 0/(Rl*T)) 

DK(5)=10**13.45*EXP(-1A./(R2*T>) 
C     DK(6)=10**14,ll*EXP(-74./(R2*T>) 

DK(7)=10**13.iS6*EXP(-89./(R2*T>> 

DK ( 8  >  * 1 0** 1 5 . 08*EXP ( -32 . / ( R2*T  > ) 
C     DK(9)=10**11.3*EXP(-51./(R2*T>> 

DK(10)--=2. 82*10**12. 0*EXP(-14.0/(R2*T>> 

DK(11)=0.0 

DK(6)=0.0 

DK(9)*0.0 
C     DK(12)=10**9.40*EXP(-8.2/(R2*T>> 
C     DK( 13) =10**7. 85*EXP(-10,0/(R2*T>> 
C     DK(14)=10**9.90*EXP(-17.1/(R2*T)) 
C      DK(15)=10**9.5*EXP(-20.8/(R2*T>) 
C     DK(li)=10**7.2*EXP(-9.2/(R2*T)) 
C     DK(17)=10**ll.l*EXP(-53.4/(R2*T>> 
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file:  toluene  Fortran  a   Kansas  state  university  um/sp  cms 


C     DK (18) =10**8. 9*EXP(-11.1/<R2*T>> 
C     DK(19)=10**9.6*EXP(-8.2/(R2*T>> 
C     DK(20)=10**11.5*EXP(-6.9/(R2*T>) 

DK  (21)  =  10**13. 0*EXP( -16. /(R2*T>) 

DK(12>=0.0 

DK(13>=0.0 

DK(14>=0.0 

DK(15>=0.0 

DK(16>=0.0 

DK(17>=0.0 

DK(18)=0.0 

DK(19)=0.0 

DK(20)-0.0 

DK(22)=10**13.0*EXP(-80.0/(R2*T>> 
C     DK(22)=0.0 

RETURN 

END 

SUBROUTINE  GRAPH1  (XMIN.XMAX,  YMIN,  YMAX,  ITITLEr  CHR) 

CHARACTER*1  CCHR,  ITITLE(76) 

DIMENSION  0X(11),0Y(11),C(51,101> 

DO  100  J=l,51 

DO  100  1=1,101 
100    C(J,I)=CHR 

XSM=XMIN 

YBG=YMAX 

XSC=(XMAX-XMIN)*0.1 

YSC=<YMAX-YMIN>*0.1 

DO  200  1=1,11 

0X(I)=XMIN+XSC*(I-1> 
200    DY(I)=YMAX-YSC*(I-1) 

XSC=10.0/XSC 

YSC=5.0/YSC 

RETURN 

ENTRY  GRAPH2  (X,Y,CHR) 

RLY=YSC*<YBG-Y>+1.5 

RLX=XSC*(X-XSM>  +  1.3 

IF  (RLX. GE. 102.0. OR. RLX.LT. 1.0)  RETURN 

IF  (RLY.0E.52.0.0R.RLY.LT.1.0)  RETURN 

C(INT(RLY),INT(RLX))=CHR 

RETURN 

ENTRY  GRAPH3 

URITE  (6,600)  ITITLE 

WRITE  (6,601)  OX 

URITE  (6.602) 

DO  300  1=1.51,5 

URITE  (A, 603)  0Y(I/5+l  > ,  (C(I,  J) ,  J=l  ,101  >  ,0Y(  1/5+1) 

Kl-I+1 

K4=I+4 
300    IF  (I.NE.51)  WRITE  (6,404)  <  (C(K,  J) , J=l , 101 ) ,K=K1 ,K4) 

WRITE  (6,602) 

WRITE  (6,601)  OX 

WRITE  (6,600) 

RETURN 

600  FORMAT  ( 1H1,25X,76A1 ) 

601  FORMAT  (/, 10X, 11 (1PE10.2) ,/) 
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file:  toluene  Fortran  a   Kansas  state  university  vm/sp  cms 

602  FORMAT  ( 14Xr '.  +  ' f20( '....  +  '),'.'  ) 

603  FORMAT  (2X. IPE10.3* '   +'il01Alr'+  'rlPE10.3) 

604  FORMAT  <14Xr ' . ' , 101A1 t ' . ' ) 
END 

SUBROUTINE  FCN<NrX. Yf YPRIME) 
DIMENSION  Y(N)fYPRIME(N> 
COMMON  DK(60) 

YPRIME<1>=-DK(1)*Y<1)-DK<2)*Y<2)*Y(1)-DK<3>*Y<2)*Y<1> 
1  +DKU0>*Y(3)*Y<7>+DK(14>*Y(8>*Y(4)+DK<15)*Y<9>*Y(5) 
1  -DK<17>*Y(1>*Y<7)-DK(21>*Y(1)*Y<8> 

YPRIME<2)=-DK(2)*Y(2)*Y(1)-DK(3>*Y(2)*Y(1)+DK(1)*Y(1) 
1  -DKC5>*Y(2>XY<4)+DK(7)*Y(4)+DK(?)*Y<6)+DK(10)*Y<3)*Y<7) 
1  -DK<11)*Y(4>*YC2>+DKU2>*Y<5>*Y(7>+DK<13>*Y<4)*YC5>+IIK(14)*Y<8>* 
1  Y  ( 4  >  +DK  <  15  >  *Y  <  9  >  *Y  <  5 )  -DK  ( 20  >  *  Y  <  1 0  >  *Y  <  2  ) 

YPRIME<3>=— DK<4)*Y<3>*Y<3)-DK<10)*Y<3>*Y<7>+DK(3>*Y(2>*Y(1> 
1  +DK<21>*Y(1)*YC8> 

YPRIME(4)=DK(2)*Y(2)*Y<1)-DK(5)*Y(2)*Y<4)-DK(6)*Y(4) 
1  -DK<7)*Y<4)-DK<11)*Y<4>*Y<2>+DK<12>*Y<5>*Y<7) 
1  -DK(13)*Y<4)*Y<5)-DK(14)*Y(8)*Y(4)-DK(16>*Y(8)*Y(4> 
I    +DK<17>*Y(1)*Y<7> 

YPRIME<S)=DK(5)*Y(2)*Y(4)-DK(7)*Y(4)-DK(18)*Y(5)*Y(?) 
I  -DK<8>*Y<5>+DK<ll>*Y(4)*YC2>-BK<12>*YC5)*Y<7>-DK'U9>*y<5>*Y<7> 
1  -DKU3)*Y<4>*Y(5>-DK<15>*Y<9>*Y<5>+DK<16>*Y<8)*Y<4) 

YPRIME  <6>=DK(8>*Y<5)-DK'<9>*YC6> 

YPRIME(7)=DK(3)*Y<2)*Y(1>+DK(S)*Y(2)*Y<4)-DK(10)*Y(3)*Y<7) 

1  +DK(11)*Y(4)*Y(2)-DK(12)*Y(5)*Y<7)-DK(17)*Y(1)*Y(7)-DK(19)*Y(5) 
1  *Y<7) 

YPRIME(8)=DK(2)*Y(2)*Y<1)-DK(14)*Y(8>*Y(4)-DK(16)*Y(S)*Y(4) 
1  -DK(21>*Y<1)*Y(8) 

YPRIME<9>=-DK(1S>*Y(9)*Y<5)+DK(16)*Y(8)*Y(4)+DK(17)*Y(1>»Y<7) 
1  +DK<19>*Y(B>*YU>+DK<21>*Y<1>*Y<8> 

YPRIME(10)=<DK'<13>*Y(4>*Y<5>-DK<20>*Y<10>*Y<2) 

YPRIME(11)=DK(4)*Y(3)*Y(3)+DK(6)*Y<4)+DK(B)*Y(5)+DK<22)*Y(1-') 

YPRIME<12>=DK<4>*Y(3>*Y<3)-DK(22)*Y(12> 

RETURN 

END 

SUBROUTINE  FCNJCNfXf YfPD> 

INTEGER  N 

REAL    Y(N)fPD(NfN)fX 

RETURN 

END 


126 


Appendix  F.   Error  Analysis 

The  calculation  of  the  energy  of  activation  (E  )  for  the  overall 
reaction  rate  as  well  as  the  product  yield  concentrations  in  the  thermal 
decomposition  of  toluene  were  performed  from  a  series  of  experimental 
measurements.   Therefore,  the  estimation  of  error  bars  for  these 
parameters  should  be  performed  using  the  propagation  of  error  technique. 

Experimental  errors  associated  with  other  variables  such  as 
temperature,  pressure  and  time  can  be  found  in  Table  F.l. 

From  the  Arrhenius  form  of  the  reaction  rate  equation,     the 
energy  of  activation  for  toluene  decomposition  is: 

EA  =  RT(S.n  A  -  in  K).  (F.l) 

where  for  first  order  reaction  K  is  given  by 

to  [T]/[T1 


K - (F.2) 

and  R  is  the  universal  gas  constant,  A  the  frequency  factor,  T  is  the 

temperature,  't'  the  reaction  time  and  [T] ,  [Tl   are  the  final  and 

o 

initial  toluene  concentrations.   The  propagation  of  error  equation  for 
the  standard  deviation  of  K  is  given  by. 


m?  °m  *  [wrf  "mo  ♦  (f)2  °l 


where  the  symbol  a  means  standard  deviation  of  term  'i'.   The  standard 
deviations  and  partial  derivatives  expressions  can  be  found  in  Tables 

F.l  and  F.2 
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Similarly,  the  standard  deviation  of  E.  is  obtained  from 

A 


2 
where  the  term  oR  from  Eq.  (F.3)  is  inserted  into  Eq.  (F.4).   In  order 

to  evaluate  Eqs.  (F.3)  and  (F.4)  the  following  set  of  conditions  was 

selected  from  the  experimental  data: 


T  =  1450  K 
[Tq]  -  2.87  k  10"7  mole/cm3 
t  =  2.0  msec 
[T]/[T]   =  0.27 

Inserting  these  values  into  the  previous  equations  a  value  of  33.06 
sec   is  obtained  for  the  deviation  of  the  constant  (corresponding  to 
5.05%  deviation)  and  a  value  of  2.5  kJ/mole  is  obtained  for  the  energy 
of  activation  which  corresponds  to  a  2.3%  deviation. 

Thus,  by  the  assumption  of  normal  distribution  95%  of  the  values  of 
E  will  be  bounded  by: 

EA  =  107  ±  4  kJ/mole  (F.5) 

The  procedure  for  the  product  yields  normalized  to  one  millisecond  is 
similar  to  the  previous  one.   In  this  case  the  equation  used  is 


Xf  Pf  1 
U  'tTp-T   «/sec)  (F.6) 

I  I 
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where  the  ratio  Xf/X  represents  the  initial  (I)  toluene  composition, 
and  the  final  (f)  product  yield,  Pf  is  the  final  pressure,  P  is  the 
initial  pressure,  and  t  is  the  reaction  time.   The  standard  deviation 
for  u  is  given  by 


_3u 

3X 


"I2 


(F.7) 


The  standard  deviation  and  partial  derivatives  for  the  different  terms 
are  found  in  Tables  F.l  and  F.2.   Using  similar  data  as  before  and 
benzene  yields  at  1450  K  a  value  of  1.13%/msec  is  obtained  for  the 
standard  deviation  which  corresponds  to  a  deviation  of  11%.   Therefore, 
the  error  bars  associated  with  the  results  for  the  product  yield  curves 
is  of  11%. 
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Table  F.l.   Standard  Deviation  of  Measured  Parameters 


Parameter 

Standard  Deviation 
a 

Comment 

t  (time) 

±0.1  msec 

estimate  from  smallest 
division  on  oscillogram 

P  (pressure) 

±10  torr 
(1.32  x  10 

estimate  from  one-half 

atm) 

of  major  division  on 

gauge. 

X,,  X 

f    0 

2.5  x  10~3 

*Calculated  from  g.c. 
analysis 

T  (temperature) 

±35 

♦Calculated  from 
reaction  time  study 

m,  [tio 

±7.18  x  10~10 

mol/cm3 

based  on  previous  value 
of  ±0.25%  for  X.  and  X 

f         0 

U2       <**/ 

2 

i      n 

n-1 


where:   X  =  variable  considered 


n  ■  number  of  replications 

0 

=  standard  deviation  of  variable  x 
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Table  F.2.   Expressions  for  Partial  Derivatives 


3K   *n  W^c 


3t        2 


3K 


3[T1      t  [T] 


3K     -1 


3[T]    t  [T] 


3  EA 

-jj-  -  R(£n  4  =  to  K) 


9EA,RT 

3K    K 


3u  =   J_^f 
3t  "  "  t2  Xj 


3u    _1 

3Xf  "  t  Xx 


3u      l_ 

9  Xf  "  ~  t  3£_ 


3u      Xf 
3  h'  "   X  2 
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Appendix  G  -  Experimental  Data 

In  this  appendix,  experimental  data  for  shock  waves  run  at  1.25% 

and  0.74%  toluene  composition  are  shown  in  Tables  Gl  through  G3.   First, 

the  data  for  shock  wave  parameters  calculated  from  shock  wave  relations 

discussed  in  Chapter  2  (Experimental  Procedure) .   These  are  shown  in 

Tables  Gl  (for  1.25  x  10"  initial  toluene  mol  fraction)  and  G2  (for 

-3 
7.4  x  10  ).   Then  Table  G.3  shows  all  the  normalized  product  yields  per 

initial  toluene  concentration  obtained  from  chromatographic  and 

computer  analysis  of  raw  data. 
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Table  G.I. 

Shock 

Datafoi 

'  toluene 

initial 
10 

mol  fraction 

of  1.25  x 

Uncorrec 

ted 

Shock 

Reaction 

Time 

T5 

P5 

[C7H8]/[C7H8]o 

(msec) 

(K) 

(atm) 

85  050501 

1.20 

1255 

9.2 

0.7586 

050901 

1.30 

1409 

11.7 

0.3346 

050902 

1.35 

1366 

11.2 

0.5130 

051001 

1.35 

1333 

10.8 

0.5186 

051301 

1.54 

1386 

11.4 

0.4226 

051302 

0.90 

1284 

11.4 

0.8808 

061101 

1.20 

1115 

10.9 

0.7443 

061201 

1.00 

1016 

9.3 

0.7856 

061301 

1.00 

1149 

11.4 

0.7839 

070601 

0.70 

977 

11.8 

0.8405 

071101 

1.00 

963 

10.4 

0.8175 

071701 

1.45 

1358 

10.7 

0.6727 

071702 

1.60 

1301 

10.7 

0.7173 

071802 

0.70 

1306 

8.8 

0.8397 

080901 

1.00 

817 

8.2 

0.8194 

100201 

1.40 

1337 

8.8 

0.6416 

100202 

1.30 

1320 

7.6 

0.7010 

100203 

1.85 

1575 

9.9 

0.1870 

100401 

1.80 

1514 

9.4 

0.2350 

100601 

1.80 

1524 

9.5 

0.1857 

100701 

1.40 

1461 

9.6 

0.3011 

100702 

1.40 

1418 

8.9 

0.3940 

100703 

1.40 

1446 

9.1 

0.2694 

100901 

1.80 

1501 

10.0 

0.2575 

100902 

1.20 

1320 

11.3 

0.7051 

100903 

0.95 

1149 

9.9 

0.7706 

101601 

1.60 

1493 

9.2 

0.1979 

101602 

1.20 

1406 

7.3 

0.4160 

101603 

1.30 

1564 

9.8 

0.2115 

101801 

1.50 

1580 

10.0 

0.1790 

102201 

1.55 

1553 

9.7 

0.1175 

102202 

1.60 

1597 

8.4 

0.1663 
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Table  G.2. 


Shock  Data  for  Toluene  Initial 
Mol  Fraction  of  7.4  x  10 


Uncorrec 

ted 

Shock 

Reaction 

Time 

T5 

P5 

[C7H8]/[C7H8]o 

(msec) 

(K) 

(atm) 

85  102801 

1.35 

1635 

9.1 

0.1418 

102802 

1.00 

1388 

9.3 

0.5134 

102803 

1.35 

1536 

9.6 

0.2467 

102901 

1.35 

1308 

9.8 

0.6343 

102902 

1.15 

1287 

10.2 

0.7883 

110401 

1.35 

1516 

9.8 

0.2381 

110402 

1.25 

1501 

10.0 

0.2434 

110403 

1.15 

1443 

10.2 

0.3728 

110501 

1.35 

1555 

7.6 

0.2421 

110702 

1.05 

1332 

11.1 

0.6247 
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Table  G.3.  Experimental  Product  Yields  Computed  by  GC's  and 

Computer  Programs. 

Corrected 

(Product  Yields/Initial  Toluene 

Composition) 

Shock 

IjClc) 

Dwell  Time 

per  msec  (%/msec) 

(msec) 

"* 

C,H, 
2  2 

C2H4 

C3H6 

Propadlene 

C3H4 

85  050501 

1255 

1.723 

0.00 

0.041 

0.076 

0.022 

0.0 

0 

050901 

1409 

1.845 

1.25 

1.759 

0.373 

0 

0.089 

0.181 

050902 

1366 

1.915 

0.00 

0.176 

0.138 

0 

0.041 

0.089 

051001 

1333 

051301 

1386 

2.100 

0.55 

1.021 

0.263 

0 

0.084 

0.151 

051302 

1284 

1.431 

0.00 

0.045 

0.081 

0.021 

0.026 

0.044 

061101 

1115 

3.525 

0.00 

0.000 

0.000 

0.005 

0.008 

0.007 

061201 

1016 

3.341 

0.00 

0.000 

0 

0 

0 

0 

061301 

1149 

3.385 

0.00 

0.003 

0 

0.004 

0.101 

0 

070601 

977 

2.830 

0.00 

0.000 

0 

0 

0 

0 

071101 

963 

3.397 

0.00 

0.000 

0 

0 

0 

0 

071701 

1358 

2.023 

0.00 

0.391 

0.123 

0 

0.045 

0.096 

071702 

1301 

2.175 

0.00 

0.255 

0.112 

0.020 

0.035 

0.066 

071802 

1306 

1.198 

0.00 

0.026 

0: 

0.006 

0.021 

0.028 

080901 

817 

3.226 

0.00 

0.000 

0 

0 

0 

0 

100201 

1337 

1.937 

0.00 

0.506 

0.222 

0.026 

0.052 

0.101 

100202 

1320 

1.-833 

0.00 

0.365 

0.228 

0.020 

0.041 

0.090 

100203 

1575 

2.451 

1.62 

7.284 

0.266 

0 

0.044 

0.106 

100401 

1514 

2.885 

3.23 

6.401 

0.476 

0 

0.071 

0.161 

100601 

1524 

2.366 

0.00 

0.000 

0 

0 

0 

0     __ 

100701 

1461 

1.971 

3.35 

7.020 

0.852 

0.044 

0.117 

0.271  " 

100702 

1418 

1.949 

1.58 

3.763 

0.799 

0.001 

0.123 

0.289 

100703 

1446 

1.946 

3.01 

7.211 

10.63 

0.030 

0.135 

0.319 

100901 

1501 

2.360 

3.16 

7.528 

0.621 

0 

0.101 

0.226 

100902 

1320 

1.751 

0.00 

0.184 

0.142 

0.019 

0.032 

0.060 

100903 

1149 

3.021 

0.00 

0.003 

0 

0.004 

0.002 

0.006 

101601 

1493 

2.147 

2.93 

9.533 

0.361 

0 

0.097 

0.190 

101602 

1406 

1.736 

2.46 

3.186 

0.972 

0.069 

0.114 

0.285 

101603 

1564 

1.853 

1.63 

9.838 

0.315 

0 

0.067 

0.160 

101801 

1580 

2.057 

2.13 

8.786 

0.288 

0 

0.068 

0.156 

102201 

1553 

2.089 

2.81 

9.755 

0.351 

0 

0.060 

0.155 

102202 

1597 

2.150 

1.56 

5.681 

0.197 

0 

0.043 

0.092 

102801 

1635 

2.102 

1.19 

6.852 

0.182 

0 

0.056 

0.166 

102802 

1388 

1.709 

0.33 

1.865 

0.484 

0.043 

0.136 

0.365 

102901 

1308 

2.101 

0.31 

0.054 

0.050 

0 

0.055 

0.151 

110401 

1516 

2.098 

2.46 

8.230 

0.447 

0 

0.088 

0.219 

110402 

1501 

1.979 

2.94 

9.442 

0.653 

0 

0.128 

0.317 

110403 

1443 

1.890 

1.71 

4.496 

0.713 

0.043 

0.129 

0.323 

110501 

1555 

2.096 

3.33 

11.695 

0.444 

0 

0.104 

0.292 

110702 

1332 

1.795 

0.00 

0.482 

0.206 

0.017 

0.062 

0.186 

135 


Corrected 

(Product  Yleld./InitUl 

.  Toluene  CcaoMltUa) 

Shock 

I  (k) 

Dwell  Time 

per  aaee 

(tjBHC) 

(msec) 

1-3  Butadiene 

C  H 
4  4 

C6H6 

Va 

Blbenxyl 

35  050501 

1255 

1.723 

0 

0 

0.52 

75.9 

0 

050901 

1409 

1.845 

0.0178 

0 

4.79 

33.5 

0.0017 

050902 

1366 

1.915 

0.0113 

0 

1.97 

51.3 

0.0140 

051001 

1333 

051301 

1386 

2.100 

0.0130 

0 

3.29 

42.3 

0.0177 

051302 

1284 

1.431 

0 

0 

0.61 

88.1 

0.0076 

061101 

1115 

3.525 

0.0039 

061201 

1016 

3.341 

0 

0 

0 

78.6 

0.0014 

061301 

1149 

3.385 

0.0010 

0.0017 

0.03 

78.4 

0.0025 

070601 

977 

2.830 

0.0012 

0 

0 

84.0 

0.0007 

071101 

963 

3.397 

0.0013 

0 

0.07 

81.8 

0 

071701 

1358 

2.023 

0.0095 

0 

2.15 

67.3 

0.0026 

071702 

1301 

2.175 

0.0072 

0 

1.78 

71.7 

0.0008 

071802 

1306 

1.198 

0.0076 

0 

0.37 

84.0 

0.0050 

080901 

817 

3.226 

0 

0 

0.03 

81.0 

0 

100201 

1337 

1.4S7 

0.0131 

0.0027 

1.83 

64.2 

0.0001 

100202 

1320 

1.-833 

0.0121 

0.0025 

1.43 

70.1 

0.0066 

100203 

1575 

2.451 

0.0061 

0 

1.12 

18.7 

0.0030 

100401 

1514 

2.885 

0.0080 

0.0010 

3.37 

23.5 

0.0031 

100601 

1524 

2.366 

0 

0 

1.48 

18.6 

0.0033 

100701 

1461 

1.971 

0.0189 

0.0043 

4.83 

30.1 

0 

100702 

1418 

1.949 

0.0317 

0.0037 

5.71 

39.4 

0.0007 

100703 

1446 

1.946 

0.0301 

0.0059 

5.17 

26.9 

0.0005 

100901 

1501 

2.360 

0.0114 

0.0025 

3.61 

25.7 

0.0010 

100902 

1320 

1.751 

0.0076 

0.0025 

0.91 

70.5 

0.0630 

100903 

1149 

3.021 

0.0016 

0 

0.03 

78.4 

0.0025 

101601 

1493 

2.147 

0.0102 

0.0034 

1.83 

19.2 

0 

101602 

1406 

1.736 

0.0349 

0.0081 

5.17 

41.6 

0.0089 

101603 

1564 

1.853 

0.0111 

0.0058 

1.19 

21.1 

0.0222 

101801 

1580 

2.057 

0.0070 

0.0034 

1.12 

17.9 

0.0048 

102201 

1553 

2.089 

0.0091 

0.0076 

1.37 

11.7 

0.0003 

102202 

1597 

2.150 

0.0061 

0 

1.41 

16.6 

0 

102801 

1635 

2.102 

0.0162 

0.0 

0.85 

14.2 

0 

102802 

1388 

1.709 

0.0673 

0.0274 

4.04 

51.3 

0 

102901 

1308 

2.101 

0.0239 

0 

1.09 

63.4 

0 

110401 

1516 

2.098 

0.0171 

0.0059 

1.61 

23.8 

0 

110402 

1501 

1.979 

0.0252 

0.0170 

2.78 

24.3 

0 

110403 

1443 

1.890 

0.0353 

0.0045 

5.16 

37.3 

0 

110501 

1555 

2.096 

0.0228 

0.0105 

2.43 

24.2 

0 

110702 

1332 

1.795 

0.0298 

0.0191 

1.89 

62.5 

0 
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Corr«cc«d 

1 

Shock 

I5(k) 

Dwell  Ilaa 

Dtphenyl 

(■«c) 

Bl phenyl 

Methane 

85  050501 

1255 

1.723 

0 

0 

050901 

1*09 

1.845 

0.0013 

0.0007 

050902 

1366 

1.913 

0.0077 

0.0071 

051001 

1333 

051301 

1386 

2.100 

0.0079 

0.0081 

051302 

1284 

1.431 

12.4  E-5 

0.0004 

061101 

1115 

3.525 

061201 

1016 

3.341 

0 

0 

! 

061301 

1149 

3.385 

0 

0.0011 

070601 

977 

2.830 

0 

0.0016 

' 

071101 

963 

3.397 

0 

0 

071701 

1358 

2.023 

0 

0.0006 

071702 

1301 

2.175 

0 

0.0003 

071802 

1306 

1.198 

6.8  £-5 

0 

080901 

817 

3.226 

0 

0 

100201 

1337 

1.437 

0 

0.0004 

100202 

1320 

1:833 

4.5  E-5 

0.0013 

100203 

1575 

2.451 

1.9  E-5 

0.0066 

100401 

1514 

2.883 

1.8  E-5 

0.0032 

100601 

1524 

2.366 

0.1  E-5 

0.0045 

100701 

1461 

1.971 

0 

0 

100702 

1418 

1.949 

0 

o 

100703 

1446 

1.946 

1.2  E-5 

0.0002 

100901 

1501 

2.360 

3.1  E-5 

0.0003 

100902 

1320 

1.751 

4.5  E-5 

0.0016 

100903 

1149 

3.021 

0.1  E-5 

0.0011 

101601 

1493 

2.147 

0 

o 

101602 

1406 

1.736 

2.9  E-5 

0.0010 

101603 

1564 

1.833 

0.0011 

0.0059 

101801 

1580 

2.057 

0.0022 

0.0042 

102201 

1553 

2.089 

0.2  E-5 

0.0004 

102202 

1597 

2.150 

0 

0- 
0 
0 
0 
0 
0 
0 
0 

0 

0 
0 
0 
0 
0 

102801 

1635   ■ 

2.102 

102802 

1388 

1.709 

102901 

1308 

2.101 

110401 

1516 

2.098 

110402 

1501 

1.979 

'  110403 

1443 

1.890 

110501 

1555 

2.096 

0 
0 
0 

110702 

1332 

1.795 
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Appendix  H.   Product  Yield-Temperature  Curves 

In  this  appendix,  the  yield  vs.  temperature  distribution  curves  for 
C_,  C,,  and  high  molecular  weight  products  are  shown.   These  are  the 
minor  products  in  toluene  pyro lysis  and  they  are:  propyne,  propadiene, 
propene,  1-3  butadiene,  vinylacetylene,  diacetylene,  biphenyl,  diphenyl 
methane,  and  bibenzyl. 
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Fig.  H. 1.   Product  Yield-Temperature  Curves  for  Minor  Products 

(a)  Propyne  Yield,  Initial  Toluene  Mol  Fraction 
(ITMF)  =  0.74  x  10"Z 

(b)  Propyne  Yield,  ITMF  =  1.25  x  10~2 

(c)  Propadiene  Yield,  ITMF  =  1.24  x  10-2 

(d)  Propadiene  Yield,  ITMF  =  0.74  x  10~2 

(e)  Propene  Yield,  ITMF  =  1.25  x  10~2 

(f)  Propene  Yield,  ITMF  =  0.74  x  10~2 

(g)  1,3  Butadiene  Yield,  ITMF  =  1.25  x  10-2 
(h)  1,3  Butadiene  Yield,  ITMF  =  0.74  x  10~2 
(i)  Diacetylene  Yield,  ITMF  =  1.25  x  10~2 
(j)  Diacetylene  Yield,  ITMF  -  0.74  x  10~2 
(k)  Vinylacetylene  Yield,  ITMF  =  1.25  x  10~2 
(1)  Vinylacetylene  Yield,  ITMF  =  1.25  x  10~2 
(m)  Biphenyl  Yield,  ITMF  =  1.25  x  10~2 

(n)  Diphenyl  Methane  Yield,  ITMF  =  1.25  x  10~2 
(o)  Bibenzyl  Yield,  ITMF  -  1.25  x  10~2 
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ABSTRACT 

A  pyrolytic  study  of  the  thermal  decomposition  of  vapor  phase 
toluene  diluted  in  argon  was  performed  by  using  the  KSU  Single  Pulse 
Shock  Tube.  The  conditions  in  this  study  included  a  temperature  range 
from  800  to  1700  K  where  only  the  portion  from  1100  to  1700  K  was 
considered  further  for  analysis  of  product  species.   The  initial  toluene 

mol  fraction  in  argon  in  the  test  section  of  the  shock  tube  was  varied 

-2  -2 

from  0.74  x  10   to  1.25  x  10   .   The  reaction  time  was  1.97  +  0.13  msec 

for  the  temperature  range  1100  -  1700  K  and  initial  toluene  mol  fraction 

-2  -1 

of  1.25  x  10   ,  2.01  ±  0.21  msec  for  7.4  x  10   initial  toluene  mol 

fraction,  and  3.25  ±  0.24  msec  for  the  temperature  range  800  -  1100  K 

-2 
and  1.25  x  10   toluene  mol  fraction.   The  reaction  pressure  (P  )  was 

varied  from  10.5  to  20  atm.   Analysis  of  the  stable  products  was 

performed  by  using  three  gas  chromatographs  containing  flame  ionization 

detectors  coupled  to  a  computer  data  station  that  reduces  the  data. 

The  experimental  results  were  divided  into  three  categories: 

1)  reaction  time  study,  2)  temperature  effect  studies,  and  3)  modeling 

of  the  experimental  data.   The  reaction  time  study  confirmed  the 

conclusions  from  other  authors  that  the  toluene  pyrolysis  is  a  first 

order  reaction.   This  conclusion  was  also  drawn  when  the  initial  toluene 

concentration  in  argon  was  varied.   The  temperature  studies  provided 

different  results  from  previous  studies  where  the  high  yield  of  benzene 

were  observed  for  the  first  time  in  low  temperature  shock  tube  studies. 

Other  species  observed  were  acetylene,  methane,  ethylene,  and  smaller 

concentration  of  propyne,  propene,  propadiene,  1,3  butadiene, 


vinylacetylene,  diacetylene,  styrene,  biphenyl,  diphenyl  methane,  and 
bibenzyl.   Other  species  present  with  high  molecular  weight  were  not 
identified  and  they  included  C„,  ^in»  C.ft  compounds.   The  modeling  of 
the  experimental  data  was  based  on  the  product  yield  observations,  and 
the  tabulated  values  for  reaction  rate  constants  obtained  from  the 
literature. 


